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ABSTRACT 


Using  the  well-known  procedures  of  Bouwkamp  (*Theoretische  en 
numerieke  behandling  van  de  bulging  door  een,ronde  opening,  Diss. 
Groningen,  Groningen  Batavia,  1941)  for  finding  the  eigenvalues  of  dif¬ 
ferential  equations  solvable  by  application  of  three-term  recursion 
formulas,  the  authors  have  developed  a  digital  computer  program  in 
symbolic  language  for  obtaining  the  eigenvalues  of  the  separated  steady- 
state  wave  equation  in  spheroidal  coordinates,  i.e.,  for  the  equation 
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where  the  parameters -f,  m,  and  hare  as  defined  by  Stratton  et  al.  in 
•^Spheroidal  Wave  Functions^ (New  York:Wiley,  1956).  The  symbolic 
language  used  is  universal  (e.g.,  add  a  to  b,  multiply  c  by  d  ,  etc.)  and 
may  be  applied  to  any  digital  computer  by  simple  conversion  to  machine 
language.  The  program  developed  requires  at  least  3400  working  cells. 
It  computes  the  eigenvalue  A,,,."  to  a  maximum  of  13  significant  places 
for  any  integer  values  of  -t  and  {-i  >  m).  The  value  of  h  is  left  arbitrary. 
The  program  also  computes  all  the  spheroidal  wave  function  expansion 
coefficients  (d  )  up  to  p  5  200, 


PROBLEM  STATUS 

This  is  an  interim  report  on  one  phase  of  the  problem;  work  is 
continuing. 
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A  DIGITAL  COMPUTER  PROGRAM  FOR  CALCULATING 
SPHEROIDAL  WAVE  FUNCTION  EXPANSION  CONSTANTS 


INTRODUCTION 

The  solution  of  the  Helmholtz  wave  equation  in  spheroidal  coordinates  leads  to 
spheroidal  wave  functions  of  the  angle  and  radial  types.  It  is  customary  to  expand  these 
wave  functions  in  two  infinite  series  consisting  of  associated  Legendre  functions  and 
spherical  Bessel  functions,  respectively.  Tables  of  expansion  coefficients  have  appeared 
in  the  literature  at  various  times.  Unfortunately,  these  tables  have  not  been  unique  since 
the  entries  depended  upon  the  definitions  of  the  functions  as  conceived  by  their  authors 
and  upon  the  author’s  choice  of  normalization.  In  the  decade  1950-1960  there  appeared, 
among  others,  two  new  books  of  tables,  one  by  Stratton,  Morse,  et  al.  (1)  and  the  other  by 
Flammer  (2).  The  first  collection  lists  the  expansion  coefficients  (hlm-t)  for  h  =  0.1 
(0.1)  1,  h  =  1  (0.2)  8,  m  =  0  (1)  8,  and  -f  =  0  (1)  8.  The  second  book  is  a  miscellany  from 
many  sources  covering  some  of  the  ranges  in  the  first  but  normalized  by  an  entirely  dif¬ 
ferent  scheme.  Flammer’s  book,  in  addition,  contains  a  cogent  account  of  wave  function 
theory.  The  extent  of  these  existing  tables  has  been  found  to  be  insufficient  for  many 
computation  needs  in  the  field  of  applied  acoustics.  The  range  of  the  eccentricity  parameter 
h  is  too  restrictive,  and  interpolation  for  h  values  not  listed  is  impractical. 

The  authors  of  the  present  report  have  therefore  turned  to  the  difficult  task  of  writing 
a  digital  computer  program  which  is  universal  enough  to  supply  spheroidal  wave  function 
expansion  constants  in  any  quantity  for  any  values  of  m,  n,  and  h.  Our  efforts  have  been 
considerably  lightened  by  the  lucid  explanation  in  Flammer’s  book,  for  which  we  acknowl¬ 
edge  a  debt  of  gratitude.  Sections  of  Flammer’s  work  appear  in  the  following  brief  explana¬ 
tion  of  the  procedure  we  used  to  formulate  the  computer  program. 


ANGLE  SPHEROIDAL  WAVE  FUNCTIONS 


Prolate  spheroidal  angle  functions  S(U,-n)  are  defined  to  be  functions  which  satisfy 
the  ordinary  differential  equation 


_d 
dr  I 


- 


S  =  0. 


(1) 


The  parameters  A  and  m  are  separation  constants.  We  specify  here  that  m  is  an  integer 
and  that  A(h)— tc-f  +  i)  as  h-o,  where  -t  is  the  order  (an  integer)  of  the  associated  Legendre 
function  For  each  value  of  the  integers  m  and  i  there  corresponds  a  function 

which  can  be  expanded  in  an  infinite  series  of  associated  Legendre  functions.  The  function 
Sn,^  is 


Sm-f  (h.-O)  =  d„  (hlm^)  (V).  (2) 

n 

If  one  substitutes  Eq.  (2)  into  Eq.  (1),  one  obtains  the  following  three-term  recursion 
formula  in  d„(hlm't); 

Pw.n  dn+2  +  %,n  ^u-2  +  d„  =  0  (3a) 


1 
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where 


=  2in  i  1)  (n  +  2m  +  2)  ^  ('ivA 

Pm,n  (2n  +  2m  +  3)  (2n  +  2m  +  5)  ’  '''  ' 

(3c) 
(3d) 

A  simple  rearrangement  of  Eq,  (3a)  gives 


n(n  -  1) 


and 


(2m  +  2n  -  1)  (2m  +  2n  -  3)  ’ 


2  (n  +  m)  (n  +  m  -  1) 


L  (2n  +  2m  +  3)  (2 


L)  -  2m^  -  r 

n  +  2m  -  1)  . 


h*  +  (n  +  m)  (n  +  m  +  1) . 


dn  ^ _ 

"^m,n  ^m,n 

If  one  multiplies  both  sides  of  Eq.  (4a)  by  P„,n-2  obtains 


N"=- 


^n-Ki^K^2 


or 


in  which 


rn.  1  +  2 


■  fa.  n-2 

n-2 


-‘nt2 
d_  ’ 


and 


ff"  =  p  „  cr 

^m,  n  — 2  m.n 

One  can  also  solve  Eq.  (4b)  for  Nn^.2  and  obtain 

- 

K  +  a],  =  — 


(4a) 


(4b) 

(4c) 


(4d) 


Equations  (4c)  and  (4d)  are  continued  fractions.  If  one  subtracts  these  equations  for  a  given 
and  the  result  is  zero,  the  number  is  a  characteristic  value  of  Eq.  (1).  In  symbols, 
to  determine  we  attempt  to  make  the  following  calculation 


[Nn+2],  -  K+aL  = 

The  first  term  on  the  left-hand  side  is  computed  by  initially  taking  n  so  large  that  is 
negligible,  and  then  finding  the  continued  fraction  [NJI+j],  for  successively  smaller  values 
of  n.  The  second  term  on  the  left-hand  side  is  computed  by  assuming 
n+2  <  0  and  then  computing  [Nn  +  j]^  for  successively  larger  values  of  n.  It  is  evident 
from  this  procedure  that  one  can  obtain  only  by  successive  approximation.  A  very 
effective  method  of  hastening  the  computation  of  A^^  has  been  devised  by  Bouwkamp  (3). 
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In  this  method  one  first  computes  A„<tby  allowing  +  1).  Then  one  computes  a 

correction  term  S(A„^)  by  taking  the  variation  of  Eq.  (4e)  due  to  a  variation  8(A„^);  in 
symbols  this  is 


Now 


K^7  +  S(A„^))  ’^n  +  2  (A„-C)  +  S(N:  +  2). 

Hence,  Eq.  (5a)  reads 

[N:  +  2  J  -  [NS  +  2  +  MN;;  +  a]j  +  »[Nn  +  2l2  =  0- 

From  Eq.  (4c)  one  finds 

(^n  +  2-Am-e  +  N-  +  ,  )2 


=  -  HK-t) 


^^2  ^  ^^*2 


Similarly,  from  Eq.  (4d) 


HNn  +  jla 


JL 


«(N™)  +  HKO 


(5b) 


(5c) 


(6a) 


=  HAu) 


1  + 


(n;;)' 


*  (n:)2  (n:_2)2  ^ 


(6b) 


Since  n  is  an  integer  it  is  convenient  to  let  n  =  't  -  m.  Substituting  Eqs.  (6a)  and  (6b)  into 
Eq.  (5c)  and  solving  for  ^(A„,.^)  one  obtains 

[^^^-™  +  2ll  -  [Ntn,+  2]2 

&(A  /)  - - 

r  /^:e-n.-2 

(6c) 


(N^-n.42)  ^  (N:t-^42)"(N^-n44)', 
'^'t-m  +  2  ^'t-m  +  2  ^•t-ni  +  4 


Repeated  corrections  by  Eqs.  (6c)  and  (5a)  lead  to  very  accurate  values  of  the  character¬ 
istic  number  A  i.  Now 

m'v 


Nn 

^ni.n  — 2 


(7a) 


The  number  N"  is  calculated  from  Eq.  (4b)  using  the  corrected  value  of  A,„.^.  The  ratio 
cln/dn_2  is  therefore  accurately  calculated.  From  Eq.  (7a)  one  deduces  that 
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=  d. 


.^n-2  ^n-4 


(7b) 


The  ejqpansion  constants  d„  can  be  normalized  in  several  different  ways.  We  adopt  here 
the  normalization  of  Morse  and  Feshbach  (4), 


(n  +  2m)!  ^  _  ('t  +  m)! 

rd  "  ('t  -  m)!  ’ 


(7c) 


Hence,  the  magnitude  of  dg  is  fixed  to  be 


do  = 


{■i  +  m)! 


f o  r  n  even .  (7d) 


(n  +  2m)! 


d„.. 


+  (2m)! 


Similarly  we  find  dj  to  be 

.  (-t+m)!  (7e) 

d,  = - - - - - for  n  odd.  '  ' 


r  '  1 

—  — 

d„  d„  ,  d, 

n  n-i  3 

d  -  d  ,  d, 

+  (1  +  2m)!> 

l^n.3 

L  n-2  n-4  IJ 

J 

All  remaining  d,,  are  found  by  substituting  Eqs.  (7d)  and  (7e)  into  Eq.  (7a). 

The  prolate  spheroidal  angle  functions  S(h,77)  can  be  converted  into  oblate  spheroidal 
angle  functions  by  making  the  transformation  h— jg.  The  prolate  expansion  coefficients 
dn(hlm^)  become  d^^c-iglm-t)  by  this  same  transformation.  Thus  to  find  d„(-jglm^)  one 
replaces  h*  inEq.  (3a),  (3b),  (3c),  and  (3d)  by  -g^.  The  normalization  procedure  is  identical 
for  both  prolate  and  oblate  cases. 


THE  COMPUTER  PROGRAM 

The  computer  program  is  divided  into  two  parts.  In  the  first  part  the  program  com¬ 
putes  the  eigenvalue  significant  places,  or  less  depending  on  the  parameters  m, 

h.  In  the  second  part  the  program  computes  any  number  of  spheroidal  wave  function 
expansion  constants  d^(h  |  m-C)  where  n  may  take  on  integer  values  up  to  n  =  400.  These 
constants  are  computed  to  13  significant  places  or  less  depending  on  the  parameters  m, 

I,  h.  The  precision  of  controls  the  precision  of  d^(h/m't). 

The  computation  of  eigenvalues  is  accomplished  by  Bouwkamp's  procedure  (3),  the 
essential  step  of  which  is  found  in  Eq.  (6c).  We  select  a  starting  value  of  A^.^  according  to 
the  given  magnitude  of  h.  A  short  discussion  of  the  method  of  selecting  a  good  starting 
value  for  A^^;  is  presented  in  Appendix  E  of  this  report.  If  this  starting  value  is  poor 
the  program  may  converge  to  an  erroneous  eigenvalue.  After  the  initial  value  of  A^  is 
set,  the  value  is  improved  by  application  of  Eq.  (6c).  Repeated  use  of  Eq.  (6c)  eventually 
leads  to  very  accurate  values  of  A^. 

When  the  computation  of  the  eigenvalue  is  completed  the  computation  of  the  constants 
is  accomplished  by  using  the  relations  found  in  Eqs.  (7a)  to  (7e). 
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The  computer  program  finally  assembled  is  reproduced  in  Appendix  B  of  this  report. 
The  meanings  of  the  various  instructions  used  in  the  program  may  be  found  in  the  publi¬ 
cation  Narec  Bulletin  No.  27  (July  27,  1961)  available  from  the  Research  Computation 
Center  at  NRL.  We  note  that  the  language  of  this  program  is  entirely  symbolic  (i.e.,  add 
a  to  b,  multiply  c  by  d ,  etc.).  The  language  is  therefore  universal.  It  may  be  converted 
to  the  machine  language  of  a  specific  computer  by  replacing  each  symbolic  command  by 
the  equivalent  machine  language  command. 

A  set  of  input  instructions  for  the  particular  version  of  this  symbolic  program  in  the 
NIP  language  of  the  Narec-NRL  is  found  in  Appendix  A. 


ADDITIONAL  REMARKS 

The  authors  of  this  report  are  in  the  process  (to  date)  of  preparing  digital  computer 
programs  for  computing  radial  and  angle  spheroidal  wave  functions  of  any  order  and  any 
argument,  using  the  eigenvalues  and  expansion  coefficients  found  by  the  program  described 
here. 
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APPENDIX  A 


INSTRUCTIONS  FOR  USE  OF  NAREC-NRL  SPHEROIDAL 
WAVE  FUNCTION  ROUTINE 


The  Narec-NRL  spheroidal  wave  function  routine  is  found  on  Narec  Tape  No.  4593 

1.  Prepare  a  machine-language  Operator  Instruction  Sheet  (see  Fig.  Al). 

2.  Prepare  a  Data  Input  Sheet  as  follows: 

a.  Insert  the  prolate  spheroidal  values  of  m,  -t,  and  h  in  successive  two-word 
entries  on  the  sheet,  reserving  two  twelve-digit  Narec  floating  point  words  (i.e.,  fractional 
part  and  exponent)  for  each  of  the  three  numbers.  These  values  will  take  up  six  lines. 

b.  Determine  the  highest  subscript  (XXX)  in  the  constants  desired.  If  l-m  is 
even,  the  subscript  will  be  even;  if  -fc-m  is  odd,  the  subscript  will  be  odd.  Insert  the  value 
of  XXX  as  a  two-line  Narec  floating  point  number  in  the  two  lines  following  the  six  noted 
above.  This  makes  eight  lines. 

c.  Count  the  actual  number  of  d’s  to  be  typed  out.  If  -t-mis  even,  the  count  is 
(XXX  +2)/2.  If  l-m  is  odd,  the  count  is  (XXX  +  l)/2.  Convert  this  count  to  a  hex  number 
and  pre-fix  zeros  to  make  up  a  “4-digit  hex  number"  (say  YYYY).  Insert  this  hex  number 
in  the  standard  word  appearing  in  the  ninth  line.  This  word  is 

YYYY  de  064a  e5. 

d.  If  l-m  is  even,  convert  XXX  +  2  into  a  four-digit  hex  number,  add  eight  zeros, 
and  insert  in  line  10. 

If  -t-m  is  odd,  convert  XXX  +  1  into  a  four-digit  hex  number,  add  eight  zeros,  and  insert 
in  line  10. 

When  this  procedure  is  complete  the  data  tape  contains  10  lines,  forming  one  block  of 
data  input.  At  the  head  of  the  block  write  “nnnn*  and  at  the  tail  write  “mminm.  • 

Any  number  of  blocks  may  be  inserted,  provided  the  separation  symbols  (nnnn  and 
mmmm)  are  used.  When  the  data  input  sheet  is  completed,  it  is  handed  to  Narec  Research 
Computation  Center  for  further  processing. 

3.  Manual  Changes 

a.  If  oblate  spheroidal  constants  are  desired,  write  on  the  blue  machine -language 
Instruction  Sheet 


Manual  changes:  line  0413- 

0587 

e7 

0587 

eb 

07b3  - 

0800 

00 

0000 

00 

07b4  - 

0001 

00 

0000 

00 
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OPEltATOri  IN37P.UCTI0N  SKKKT 
RCC  FYoblem  Number 

Problem  Title  _ 

Run 


Date  

Estimated  Computation  Time  /S  /i/AJ 

Programmer 

Tel.  Ext.  _ 


OUTPUT  FQF^J^AT:  Copies  original  carbons. 

Ar  He»«deoiiRiU.j  Decimalt  ^  cols.  57  Rexadecimal,  Decimal:  —  cols. 


INPU' 


TAPE  INFORMATION 


Tape 

No. 

Locations 
From  Thru 

Check  Sum 
( if  known) 

L 

Tape 

No. 

Locations 
From  Thru 

Check  Sum 
(if  known) 

1 

04‘0 

IB 

2 

OCoo 

B51 

IB 

3 

El 

^  -M 

B 

L 

B 

5. 

m 

SPECm  I^JPUT  INSTRUCTIONS.  MANUAL  CHANGES.  AND  STARTING  ORDER 


f)  /  ,-•  tAj 


c  ' 

y<>. 


3  V/3 

o7^i 


0  68  7  C7  os-.7e7- 

O  ?  Co  oC 


67  ■  f 


O  !  0  do  0  a  O  Oo  0 


PRIHT-0UT3 


starting  Order  LO  RO  _ 

SPECIAL  IHS'inUCHOHS,  ADDITIOmL  ItJFORmTION.  ETC. 


Fig.  A1  -  Sample  of  Operator  Instruction  Sheet 
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b.  K  the  tolerance  on  is  to  be  changed,  compute  the  new  tolerance  as  a  two- 
lifte  Narec  floating  point  number  and  write  this  number  as 

Manual  Change;  line  0610 - 

0611 _ 

The  tolerance  on  is  fixed  at  IQ-n  in  the  program. 

c.  If  the  tolerance  used  in  comparing  '''**^*^  "^n  +  to  be  changed, 

compute  a  new  tolerance  as  a  two-line  Narec  floating  point  number  and  write 


Manual  Change:  line  OSfe; 

05f  f  : 


This  tolerance  is  started  at  10“^  in  the  program. 

d.  If  the  eigenvalue  tolerance  (=  I  tol)  comparing  the  change  computed  eigenvalue 
for  increasing  N’s  is  to  be  changed,  compute  a  new  tolerance  as  a  two-line  Narec  floating 
point  number  and  write 


Manual  Change;  line  06(10: 

06dl; 


This  tolerance  is  fixed  at  10"'*  in  the  program. 

e.  If  the  nth  ratio  tolerance  (=  n  tol)  comparing  the  change  in  n!?_„+4  as  the 
number  of  n’s  increase  is  to  be  changed,  compute  a  new  tolerance  as  a  two-line  Narec 
floating  point  number  and  write 

Manua  1  change ;  line  06C3: - - 

06C4:  _ 


This  tolerance  is  fixed  at  10"''  in  the  program. 

f.  Starting  values  for  depend  upon  whether  h  is  considered  “small”  or  “large.” 
A  short  discussion  of  the  method  of  making  a  proper  selection  of  the  starting  value  for 
A^^  is  found  in  Appendix  E. 

As  an  example  of  the  use  of  the  instructions  outlined  above,  we  select  the  following 
as  an  appropriate  example  of  the  preparation  of  a  data  input  sheet. 

Example:  Let  m  =  1,  -t  =  7,  and  h  =  10  for  the  oblate  case.  We  desire  the  largest 
expansion  constant  to  be  (note  that  the  subscript  must  be  an  even  integer).  The  data 
input  sheet  reads  as  follows. 


NAVAL  RESEARCH  LABORATORY 


9 


n  n  n  n  n 

080000  000000 
000100  000000 
070000  000000 
000400  000000 
OAOOOO  000000 
000400  000000 
OfOOOO  000000 
000600  000000 
001 fde  064Ae5 
O03EOO  000000 

m  m  m  m  m 

This  sheet  is  sent  to  Narec  Research  Computation  Center  for  further  processing. 

Since  the  spheroid  is  oblate,  make  the  following  “manual  changes*  on  the  standard 
blue -colored  Narec  Operator  Instruction  Sheet: 

line  0413-  0587e7  0587eb 

07b3  -  080000  000000 


07b4  -  000100  000000 


APPENDIX  B 


REPRODUCTION  OF  COMPUTER  PROGRAM  IN  SYMBOLIC  FORM 


NOTE:  Any  group  of  letters  and  numbers  not  specifically  appearing  in  the  list  of  commands 
of  Narec  Bulletin  #27  is  to  be  considered  a  machine  address.  For  example,  spro,  kl,  k2, 
etc.  are  addresses,  while  p,  pd,  mp  are  standard  commands.  \  Method  of  reading:  at 
address  spro  in  the  computer  memory  place  the  command  “go  to  left  half  of  address  bgn 
etc. 


spro  lo  bgn 
0 

kl  kl  a 

ro  512 

k2  J3  p 

-J2  pd 
P  J9 
J3  p 
J5  pd 
P  Jio 
J3  P 
mp  j6 
P  Jll 
-J5  pd 
P  J12 

J8  pd 
P  J13 
j2  p 
mp  J2 
mp  J6 
P  J14 
mp  J6 
P  J15 
J2  p 
J5  pd 
P  J16 
J2  p 
mp  j6 
P  J17 
-j5  pd 
P  J18 
j6  pd 

P  J19 

-j8  pd 
P  j20 
J8  pd 
j6  pd 
P  J21 
J4  p 
mp  J4 
P  J22 
P/  J6 
P  j23 


J22  p 
mp  j22 
P  J24 
J3  p 
mp  jlO 
P  j25 
Jll  p 
mp  JIO 
-jl4  pd 
-J5  pd 
P  j26 
jl2  p 
mp  J13 
P  J27 
J26  p 
P/  J27 
mp  J22 
j25  pd 
kfdg  pd 
P  J28 
fl  yadl 
k3  J9  p 

p  JlOO 
k4  JlOO  p 

cmp  J5 

sn7  fl  k5 

fl  k6 
sn8  fl  k7 

k7  JlOO  p 

-J6  pd 
p  JlOO 
fl  k4 

k5  jq20  p 

P  Jq30 
P  Jq3i 
fl  kl2 
k6  J5  p 

P  Jq30 
p  Jq3i 
fl  kl2 
k8  J2  p 

Jq31  pd 
P  J32 
J5  pd 
mp  J32 
P  J33 
J15  p 
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-j5  pd 
P  j34 
jq31  P 
mp  j6 
P  j35 
jl8  p 
j35  pd 
P  j36 
j8  pd 
mp  j36 
P  337 
334  p 
p/  337 
P  338 
35  P 
-j38  pd 
mp  323 
333  pd 
P  339 

k9  fr  0 

klO  3q31  p 

-35  pd 
p  340 
317  p 
3q31  pd 
P  341 

-35  pd 

P  342 
3q31  pd 
P  345 

36  pd 

P  344 

-38  pd 

P  343 
3q3i  p 
mp  340 
mp  341 
mp  342 
p/  345 
p/  345 

P/  343 

p/  344 
mp  324 
P  3101 
kll  fr  0 

kl2  35  P 


P  3900 

39  P 

38  pd 

khta 

P  3192 
3q30  p 
p  3103 

P  3q3i 

kl3 

3103  p 

36  pd 

snl 

p  3193 

cmp  3192 

sn9 

fl  kgl4 
fl  kl9 

sna 

fl  kl9 

kgl4 

nl  kl4 

kl4 

kl6  a 
lo  kl5 

kl5 

lal  kll 
il  kll 

kl  6 

kl6  a 
lo  512 

kl7 

fl  klO 

3101  p 

kl8 

p  3490 

adl  2 
kl8 

3q3i  P 

36  pd 

P  3q3i 
fl  kl3 

kl9 

3q30  p 
p  3q3i 
p  3104 

kgl9 

3104  p 

36  pd 

sn2 

p  3104 
cmp  3102 

snb 

fl  qql 
fl  kg24 

snc 

fl  kg24 

qql 

nl  k20 

k20 

k22  a 
lo  k21 

k21 

lal  k9 
il  k9 

k22 

k22  a 
lo  512 
fl  k8 

339  P 

k23 

p  3300 

adl  2 
k23 

3q3i  p 

36  pd 
p  3q3i 
fl  kgl9 

kg24 

3q30  p 
p  3q3i 
p  3195 

k24 

3105  p 

36  pd 

sn3 

p  3105 

cmp  3102 

snd 

fl  k25 

gk24 

fl  hvlO 
fl  hvlO 

k25 

-3400  p 

k26 

p/  3900 

k27 

-3300  pd 

328  pd 

k28 

P  3901 

P  3IO6 
adl  2 
k25 
adl  2 
k26 
adl  2 
k27 
adl  2 
k28 
fl  k24 

hvlO 

39  P 

36  pd 

p  3120 
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hvl  1 

P  J123 
-j6  pd 

P  jl24 
Jq3C  p 

P  jl21 
jl21  p 

sn4 

j6  pd 

P  jl21 

sne 

cmp  jl23 
fl  hvL2 

snf 

fl  hvl5 
fl  hvl5 

hvl2 

j401  p 

hvl3 

p/  j901 

wblO 

P/  j901 

hvl4 

p  vlOO 

hvl5 

adl  2 
hvl2 
adl  2 
hvl3 
adl  2 
hvl4 
adl  2 
wblO 
fl  hvll 
Isl  hvl4 

hvl7 

hvl8 
jS  p 

hvia 

P  0 

hv22 

vlOO  p 

hw21 

cmp  j5 

sng 

fl  hw22 
fl  hv28 

hw22 

fl  hw22 
Isl  hv22 

hu22 

hv23 

Isl  hv22 

hw23 

hw24 
adl  2 

hw24 

hv23 

0  p 

hw25 

p  wlOO 
jq30  p 

P  jl22 
wlOO  p 

hv23 

mp  0 

hv24 

p  wlOO 

hv25 

jl22  p 

sn5 

j6  pd 
p  jl22 

snh 

cmp  jl24 
fl  hv26 

sni 

fl  hv27 
fl  hv27 

hv26 

adl  2 

hv27 

hv23 
fl  hw25 
adl  2 

hv22 
adl  2 
hv24 
adr  2 
hw24 
adl  2 
hw25 

jl24 

jl24  p 
-j6  pd 

P  jl24 
fl  hv22 

0 

hv28 

0 

0 

0 

Isr  hv24 

hv29 

hv29 
j5  p 

hv30 

P  0 

WlOO  p 

snj 

cmp  j5 
fl  hv31 

snk 

fl  hv33 
fl  hv31 

hv31 

yzl2  p 

hv32 

WlOO  pd 

hv33 

p  yzl2 
adl  2 
hv30 
adl  2 
hv3  2 
fl  hv30 
yzl2  p 

k29 

j5  pd 
p  yzl2 
j9  P 

k30 

P  jl07 
jl07  p 

k31 

j6  pd 

P  jl07 

P  Jill 

P  jll2 
p  jq3i 
nl  k31 
k33  a 

k32 

lo  k32 
lal  kll 

k33 

il  kll 
k33  a 

k34 

lo  512 
fl  klO 

k35 

nl  k35 
k37  a 

k36 

lo  k36 
lal  k9 

k37 

il  k9 
k37  a 

k38 

lo  512 
fl  k8 

sn6 

J39  P 
-j28  pd 
p  jl08 
mp  wzlO 
/P/ 

p  jl09 
-jlOl  p 
p/  jl08 

snm 

/P/ 

cmp  jl09 

snn 

fl  kg39 
fl  kg39 

fl  k30 
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kg39 

jq30  p 
p  jq3i 

k39 

j9  p 
-j6  pd 

P  jiio 

k40 

Jill  P 
-jS  pd 

P  jl  11 

k41 

jllO  p 
-jlll  pd 
fcl  k55 

k42 

jll2  p 
p  jq3i 
nl  k43 

k43 

k45  a 
lo  k44 

k44 

lal  kll 
11  kll 

k45 

k45  a 
lo  512 

k46 

fl  klO 
nl  k47 

k47 

k49  a 
lo  k48 

k48 

lal  k9 

11  k9 

k49 

k49  a 
lo  512 

kSO 

fl  k8 
j39  p 
-j28  pd 
jll3  pd 

P  jll4 
-jlOl  p 
p/  jll4 

P  J113 

k51 

jll3  p 
mp  jli3 
P/  jioi 

k52 

p  k300 

k53 

jll2  p 
-j6  pd 

P  jll2 

k54 

adl  2 
k52 

fl  k40 

k55 

Isr  k52 
k56 

k56 

j5  p 

P  0 

hv56 

j5  P 
p  yzlO 

hv57 

k300  p 
cmp  j5 

sno 

fl  k58 
fl  hvG4 

snp 

fl  k5S 

k58 

k300  p 
cmp  j5 

snq 

fl  k59 
fl  kv63 

snr 

fl  k59 

k59 

yzlO  p 

k60 

mp  k300 

k61 

p  ylOO 

p  yzlO 
adl  2 
k58 
adl  2 
k60 

fl  k58 

kv63  adl  2 
hv57 

Isl  hv57 
k58 

Isl  hv57 
k6G 
adl  2 
k61 

fl  hv56 
hv64  Isr  k61 
kv65 
kv65  j5  p 

P  0 

k64  ylCO  p 

cmp  j5 
sns  fl  k65 

fl  k67 
snt  fl  k65 

k65  yzll  p 

k66  ylOO  pd 

p  yzll 
adl  2 
k64 
adl  2 
k66 

fl  k64 
k67  yzl2  p 

yzll  pd 
p  yzl3 
jU3  p 
-jlC6  pd 
P''  yzl3 
p  yzl4 
j28  pd 

snu  p  j28 

yzl4  p 

/p/ 

-jl26  pd 
k62  fcl  k79 

fl  ltd 
k70  pt  11 

f  1  hov  4 
k68  TO  k68 

lo  dj8 
k69  jl51  p 

p  jl03 
P  jl04 
P  jl05 
p  jl21 
P  jl22 
P  jl23 
p  yzll 
P  yzl2 
P  J107 
P  jlll 
P  jll2 
P  jll3 
nl  apne 
dj8  dj8  a 
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ro  512 
dj9  j5  p 

P  j900 
djlO  dj50  p 

Isl  djlO 
gk24 
j9  P 
j8  pd 
p  jlU2 
jl  51  p 
p  jl03 
p  jl04 
p  jl55 
Isl  wl8 
kl8 

Isl  w23 
k23 

Isl  wl8 
k25 
Isl  w9 
k26 

Isl  w23 
k27 

Isl  w91 
k28 

fr  khb 

djll  Isr  djf6 
djl2 
djl2  j5  p 

p  0 

jq30  p 
p  wylO 
vhlO  wylO  p 
j6  pd 
p  wylO 
vhll  j901  p 
'P/ 

tiov6  cmp  j5 
fl  vhl2 
fl  vhl6 
fl  vhl2 
vhl2  wylO  p 
mp  j6 
jl7  pd 
j5  pd 
p  wyll 
-j6  pd 
mp  wyll 
p  wyl2 
wylO  p 
jl7  pd 
p  wyl3 
-j5  pd 
mp  wyl3 
p  wyl4 
wyl2  p 
p/  wyl4 
P/  j22 

vhl3  mn  jOOl 
vhl4  p  blOO 
vhl5  adl  2 
vhll 
adl  2 
vhl3 
adl  2 


vhl4 
fl  vhlO 

vhl6 

Isr  vhl4 
vhl7 

vhl7 

j5  p 

P  0 

vhl9 

cxl  0  p 
j5  pd 
p  cxl  0 

vh20 

blOO  p 

/p/ 

cmp  j5 
fl  vh21 
fl  vh24 
fl  vh21 

vh21 

cxlC  p 

vh22 

mp  blOO 

vh23 

p  clOO 
p  cxlO 
adl  2 
vh20 
adl  2 
vh22 
adl  2 
vh23 
fl  vh20 

vh24 

Isr  vh23 
vh25 

vh25 

j5  P 

P  0 

vh27 

jq30  p 
p  cqlU 

vh28 

cqlO  p 
j6  pd 
p  cqlO 
cqlO  p 
jl7  pd 
p  qclO 
p  qcl2 

vh29 

clOO  p 

/p'' 

cmp  j5 
fl  vh30 
fl  vh35 
fl  vh30 

vh3u 

cqlO  p 
/P/ 

cmp  qclO 
fl  vh31 
fl  vh32 
fl  vh32 

vh31 

qclO  p 
-j5  pd 
p  qcll 
mp  qcl2 
p  qcl2 
qcl  1  p 
p  qclO 
fl  vh30 

vh32 

qcl2  p 
p/cqlO 

vh33 

mp  ClOO 

vh34 

p  vxlO 
adl  2 
vh29 
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adl  2 
vh33 
adl  2 
vh34 
fl  vh28 

vh35  Isi'  vh34 
vh36 
vh36  j5  p 
P  0 

vh38  vxlO  p 
/p/ 

cmp  jS 
fl  vh39 
fl  vh80 
fl  vh39 
vh39  qq2  p 
vh40  vxlO  pd 
p  qq2 
adl  2 
vh38 
adl  2 
vh40 

jhb  fl  vh38 
j9  P 

cmp  jl51 
fr  ersp 
fr  vh44 
p  bclO 
vh41  -j5  pd 
p  bell 
cmp  jl51 
fl  vh42 
fl  vh42 
bclO  p 
mp  bell 
p  bclO 
bcl  1  p 
fl  vh41 
vh42  j3  p 
j2  pd 
Cmp  jl51 
fr  ersp 
fl  vh54 
vh43  p  bcl2 
-j5  pd 
p  be  13 
cmp  jl51 
fl  vh55 
fl  vh55 
bcl  2  p 
mp  bcl3 
p  bcl2 
bcl3  p 

vh44  fr  vh43 
j5  P 
p  bclO 
fl  vh42 
vh54  j5  p 

p  bcl2 
vh55  be 12  p 

p/bclO 
P/  qq2 
p  qq3 
p  vvxl 
vh45  clOO  p 


/P/ 

cmp  j5 
fl  vh46 
fl  vh48 
fl  vh46 

vh46 

clOO  p 
mp  qq3 

vh47 

p  vvx2 
adl  2 
vh45 
adl  2 
vh46 
adl  2 
vh47 
fl  vh45 

vh48 

rep  6 
nr  eop 

vh49 

0 

0 

vh50 

vvxl  p 

vh51 

pt  11 
adl  2 
vh50 
fr  vh49 

j2 

080000 

0 

000100 

0 

j3 

080000 

0 

000200 

0 

j4 

080000 

0 

000100 

0 

jl50 

OcOOOO 

0 

4  zz 

0 

j5 

O8C000 

0 

1  zz 

0 

j6 

080000 

0 

2  zz 

0 

J7 

OcOOOO 

0 

2  zz 

0 

j8 

080000 

0 

3  zz 

0 

j9 

0 

0 

0 

0 

jio 

0 

0 

0 

0 

jll 

0 
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0 

0 

0 

jl2  0 

0 
0 
0 

jl3  0 

0 
0 
0 

Jl4  0 

0 
0 
0 

jl5  0 

0 
0 
0 

jl6  0 

0 
0 
0 

jl7  0 

0 
0 
0 

jl8  0 

0 
0 
0 

J19  0 

0 
0 
0 

J20  0 

0 
0 
0 

j21  0 

0 
0 
0 

j22  0 

0 
0 
0 

j23  0 

0 
0 
0 

j24  0 

0 
0 
0 

j25  0 

0 
0 
0 

j26  0 

U 
0 
0 
0 
0 


0 

0 

j28  0 

0 
0 
0 

jlOO  0 
0 
0 
0 

jq20  0 
0 
C 
0 

jq30  0 
0 
0 
0 

jq31  0 
0 
0 
0 

j32  0 

0 
0 
0 

j33  0 

0 
0 
0 

j34  0 

0 

0 

0 

j35  0 

0 
0 
0 

j36  0 

0 
0 
0 

j37  0 

0 
0 
0 

j38  0 

0 
0 
0 

j39  0 

0 

0 

j40 


J41 


j42 


J27 


ocoooooooooo 
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0 

j45  0 

0 
0 
0 

j44  0 

0 
0 
0 
0 

j43  0 

U 
0 
0 

jlOl  c 

0 
0 
0 

jl02  0 
0 
0 
0 

jl03  G 

0 
0 
C 

jl04  0 

0 
0 
c 

jl05  0 

0 
0 
0 

jl06  0 

0 
0 
0 

jl2G  C 


jl21  0 
0 
0 
0 

jl22  0 
0 
U 
0 

jl23  0 
0 


yzl2  0 
0 
0 
0 

jl07  0 
0 
0 
0 

Jill  0 
0 
0 


0 

jll2  0 
0 
0 
0 

jl08  0 
0 
0 
0 

WzlO  0a3d70 
a3cl70a 
f  f  f  aOO 

jl09  G 
0 
0 
0 

jlio  0 
0 
0 
0 

jll3  0 
0 
0 
0 

jll4  0 
0 
0 
0 

yzlO  0 
0 
0 
0 

yzll  0 
0 
0 
0 

yzl3  0 
0 
0 
0 

yzl4  0 
0 
0 
0 

jl26  OafecO 
0 

ffdcOO 

0 

wylO  0 
0 
0 
0 

wyll  0 
0 
0 
0 

wyl2  0 
0 
0 
0 

wyl3  0 
0 
0 
0 
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wyl4  0 
0 
0 
0 

cxlO  0 
0 
0 
0 

cqlO  0 
0 
0 
0 

qclO  0 
0 
0 
0 

qcl2  0 
0 
0 
0 

qcll  0 
0 
0 
0 

qq2  0 
0 
0 
C 

bcl3  0 
0 
0 
0 

bclO  0 
0 
0 
0 

bcl2  0 
0 
0 
0 

bcl  1  0 

0 
0 
0 

qq3  0 
0 
0 
0 

jl51  0 

0 
0 
0 

wl8  j400  zz 

0 

w23  j300  zz 

0 

w9  j900  zz 

0 

w91  j901  zz 

0 

bgn  rcl  psm 

lo  i4U0 

400  t  psm 

j2  8 
t  prep 


vh48  1 
040083 
240030 
lo  kl 

psm  0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

prep  0 
0 

wj71  0 
0 

apne  240081 
0400fd 

apnf  apnf  a 
ro  512 
jq30  p 
P  jq31 

eop  fl  kl2 
stop 

nrun  240083 
040030 
lo  bgn 
0 

fset  pc  37 
pc  37 
pc  7 
pc  4 
pc  4 
j2  p 
0300fl 
pc  37 
pc  6 
pc  4 
pc  4 
j3  p 
0300f 1 
pc  37 
pc  8 
pc  4 
pc  4 
j4  p 
0303fl 
pc  37 
pc  37 
pc  37 
pc  16 
pc  16 
pc  16 
pc  16 
pc  41 
j28  p 
fl  k70 
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vh80  J17  p 

cmp  jl51 
fr  ersp 
fr  vh83 
vh81  fl  vh84 
-j5  pd 
p  hall 
cmp  jl51 
fr  vh82 
fr  vh82 
halo  p 
mp  hall 
p  ha  10 
hal  'i  p 

vh82  fr  vh81 
qq2  p 
ha 10  pd 
p  qq2 

ersp  fr  jhb 
nl  vh83 
vh83  0bad82 
j5  p 
p  halo 
fr  vh82 
halo  0 
0 
0 
0 

hall  0 
0 
0 
0 

vh84  jl7  p 

jq30  pd 
p  halO 
fr  vh81 
djSO  jq30  p 
P  jq3i 
jl50  p 
P  Jill 
dj42  p  jll2 
jl51  p 
p  jll3 
j5  P 

dj41  p  ggll 

dJ51  j9  p 

p  jllO 
dj52  jlll  p 
-j6  pd 
P  Jill 
dj53  jllO  p 

-jlll  pd 
d j40  fcl  d j74 
dj54  jll2  p 
P  jq3i 
nl  dj55 
dj55  dj57  a 

10  dj56 
dj56  lal  kll 

11  kll 
dj57  dj57  a 

lo  512 
dj58  fl  klO 
nl  dj59 
dj59  dj61  a 


10  dj60 
dj60  lal  kO 

11  k9 

dj61  dj61  a 
lo  512 

dj62  fl  k8 
j39  P 
-j28  pd 
jll3  pd 
P  jll4 
-jlOl  p 
p/jll4 
P  jll3 
dj63  p  gglO 
p  nlst 
adl  4098 
dj63 
jll2  p 
-j6  pd 
P  jll2 
fl  dj52 

dj64  Isl  wj65 
dj40 
jll2  p 
j6  pd 
P  jlll 
jq30  p 
P  jiio 
fl  dj52 
wj65  dj65  zz 
0 

dj65  adl  2 
dJ63 

Isl  dj63 
dj70 

dj70  0  p 

cmp  j5 
fl  djf6 
fl  dj72 
djf6  p  j901 
skip 
adl  2 
djf0 
adl  2 
dj70 
fl  dj70 
0 

dj72  Isl  wjn 
dj71 
Isr  w9 
dj73 

dj71  adr  0 
dj73 
dj73  j5  p 
P  0 

fl  djll 
k79  jl60  p 
j5  pd 

hov7  p  J160 
cmp  j  7 
fl  k69 
fl  k69 
fl  kv70 
jl60  0 
0 
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kv70 


kten 


dJ74 


wj74 

wj41 

wJ63 

wise 

nlst 

nlsu 

cij75 


0 

0 

wzlO  p 
p/kten 
p  wzlO 
jl51  p 
p  jl60 
f IhovS 
080000 
0 

2  zz 
0 
0 
0 
0 
0 

nlst  p 

/P/ 

p  su 
Isl  wj74 
d  j40 

Isl  wj41 
dj41 

Isl  wj63 
dj63 
jl50  p 
kten  pd 
p  wl50 
P  Jill 
fl  dJ42 
dJ75  zz 
0 

ggll  zz 
0 

gglO  zz 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

-nlsu  p 
/nlst/pd 

/p/ 

emp  ntol 
fl  dj64 
fl  dj64 
nlst  p 

/P/ 

p  nlsu 
Isl  wJ41 
dJ41 

Isl  wj63 
dj63 
wl50  p 
kten  pd 
p  wl50 
P  jlll 


fl  dj42 
ntol  OafecO 
0 

ffdcOO 

0 

ltd  j28  p 
/P/ 

P  j29 
Isr  mbl 
k62 

fl  k79 

mbl  ltc2  zz 
0 

ltc2  j28  p 

/P/ 

-j29  pd 
/P/ 

-Itol  pd 
fcl  ilOO 
fl  fset 
0 

100  j28  p 

/p/ 

P  j29 
fl  k79 
Itol  OafecO 
0 

ffdcOO 

0 

j29  0 

0 
0 
0 

kfdg  080000 
OOOOfO 
1  zz 
0 

yadl  nrlb  p 
-j5  pd 
fcl  zin 
fl  vav 
vav  j4  p 

/P/ 

emp  m992 
fl  k3 
fl  gnul 
fl  gnul 
zin  j4  p 

/P/ 

emp  m992 
fl  k3 
f 1  gnu5 
fl  gnu5 
gnul  j3  p 

-j2  pd 
mp  j6 
j5  pd 
p  m9r2 
mp  m9r2 
p  m9r3 
mp  m9r2 
p  m9r4 
mp  m9r2 
p  m9r5 
mp  m9r2 
p  m9r6 
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mp  m9r2 
p  mSr7 
mp  m9r2 
p  m9r8 
j4  p 
mp  j4 
p  m9ra 
mp  j4 
p  m9rb 
mp  j4 
p  m9rc 
mp  J4 
p  m9rd 
j2  p 
mp  j2 
p  m9re 
mp  mOre 
p  mare 
gnu2  j4  p 

mp  m9r2 
p  m9rf 
m9r  3  p 
mSrO  pd 
p/  m8rl 
p  m9rg 
m8re  p 
mp  m8r3 
p  m3rh 
-mCrh  p 
m9r3  pd 
m8r2  pd 
mp  m9r2 
p/  mSr4 
P/  j4 
p  m9rj 
m9r3  p 
mp  mSrS 
m9r5  pd 
mSr6  pd 
mp  mSrO 
p  m9rk 
m9r3  p 
j5  pd 
mp  m9re 
mp  m8r8 
p  m8rh 
-m9rh  p 
m9rk  pd 
p/  m8r7 
p/  mCra 
p  m9rn 
m9r6  p 
mp  mSra 
p  m9rh 
m9r4  p 
mp  m8rb 
p  m9rk 
m9r2  p 
mp  m8rc 
m9rk  pd 
m9rh  pd 
p/  m8r9 
p/  mSr9 
p  m9rq 
ra9r2  p 


mp  m8re 
p  mCrh 
mDr4  p 
mp  mSrd 
m9rh  pd 
mp  m9re 
p/  mSrG 
p  m9rt 
m9r2  p 
mp  mare 
p/  mSrl 
-ra9rt  pd 
m9rq  pd 
p/  mGrb 
p  mCru 
m9r7  p 
mp  m8rg 
p  mGrh 
m9r5  p 
mp  m88f 
p  m9rk 
m9r3  p 
mp  clal 
m8rh  pd 
m9rk  pd 
m9rh  pd 
p/  mSrf 
p/  mSrf 
p  m9rt 
m9r5  p 
mp  mSrn 
p  mGrh 
mGr3  p 
mp  mSrj 
p  mGrk 
m8rp  pd 
m9rh  pd 
mp  mGre 
p/  m8rk 
p  m8rq 
m9r3  p 
j5  pd 
mp  mare 
mp  m991 
p/  m8rl 
-m9rq  pd 
m9rt  pd 
p/  m9rc 
p  m9rw 
m9r8  p 
mp  m8rq 
p  mGrh 
m9r6  p 
rap  m8rs 
p  m9rk 
m9r4  p 
mp  m8rt 
p  m9rq 
m9r2  p 
rap  m8ru 
m9rq  pd 
m9rk  pd 
m9rh  pd 
p/  m8r7 
p/  m8r7 
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p  m9rh 
mOrG  p 
mp  mSrv 
p  mOrk 
m9r4  p 
mp  mGrw 
p  mOrq 
mSr2  p 
mp  mSry 
mOrq  pd 
mSrk  pd 
mp  mLre 
p/  m8r3 
p/  m8r7 
p  m9rk 
m8r  z  p 
mp  m9r4 
p  m9rq 
m9r2  p 
mp  mOrO 
m9rq  pd 
mp  mare 
p  f  m8rk 
p  mOrq 
m9r2  p 
mp  m9re 
p/  m9rl 
p  mSrt 
-m9rt  p 
m9rq  pd 
-m9rk  pd 
m9rh  pd 
p/  m9rd 
p  mOrx 
f  1  hov  1 
gnu 3  mOrf  p 
m9re  pd 
-m9rg  pd 
-m9rj  pd 
-mOrn  pd 
-m9ru  pd 
-mOrw  pd 
-m9rx  pd 
gnu4  p  j28 
fl  k3 

gnu5  j4  p 
mp  j4 
P  g3f0 
mp  j4 
p  g3fl 
mp  J4 
p  g3f2 
mp  j4 
P  g3f3 
j3  P 
-j2  pd 
mp  JG 
j2  pd 
j5  pd 
P  g3f4 
mp  g3f4 
P  g3f5 
mp  g3f5 
P  g3f6 
mp  g3f5 


P  g3f7 
j2  p 
mp  j2 
P  g3fS 
mp  g3f8 
P  g3fG 
mp  g3f8 
P  g4f0 
gnuG  g3f4  p 
mp  j4 
mp  jG 
p  g4fl 
g3f5  p 
-g3fS  pd 
j5  pd 
P/  jS 
P  g4f2 
p/  g2ck 
P/  j4 
mp  g3f4 
P  g4f3 
g3f6  p 
mp  m8r0 
p  gSaO 
g3f5  p 
mp  m9D2 
P  gSal 
g3f5  p 
mp  m991 
j5  pd 
rap  g3f8 
mp  jG 
p  g5a2 
-g5a2  p 
g3f9  pd 
J5  pd 
g5al  pd 
gSaO  pd 
p/  m8r4 
p/  g3fO 
P  g4f4 
g3f6  p 
rap  ra8ra 
p  gSaO 
g3f5  p 
mp  g2c2 
p  gSal 
g3r5  p 
mp  g2cl 
g2c4  pd 
mp  g3f8 
mp  j6 
P  g5a2 
g3f9  p 
mp  g2c5 
-g5a2  pd 
m8rd  pd 
gSal  pd 
gSaO  pd 
rap  g3f4 
p/  m8rk 
P/  g3fl 
P  g4f5 
g3f7  p 
mp  m8rg 
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p  g5a0 
B3fG  p 
mp  e2c3 
p  g5al 
g3f5  p 
mp  g2c6 
p  g5a2 
g3r6  p 
mp  m092 
p  g5a3 
g3f5  p 
mp  g2c 1 
m991  pd 
g5a3  pd 
mp  g3f8 
mp  m992 
p  g5a3 
g3f5  p 
mp  g2c5 
g2c8  pd 
mp  g3f9 
mp  m99 ! 
p  g5a4 
g4f0  p 
mp  j6 
p  g5a5 
-gSaS  p 
g5a4  pd 
-g5a3  pd 
g2c7  pd 
g5a2  pd 
gSal  pd 
gSaO  pd 
p/  mSr? 
p/  g3f2 
p  g4f3 
g3f7  p 
mp  m8rq 
p  gSaO 
g3f6  p 
mp  g2c9 
p  g5al 
g3f5  p 
mp  g2ca 
p  g5a2 
g3fG  p 
mp  g2cc 
p  g5a3 
g3f5  p 
mp  g2cd 
g2ce  pd 
g5a3  pd 
mp  g3fC 
P  g5a3 
g3f5  p 
mp  g2cf 
g7cg  nd 
mp  g3f9 
n  g5a4 
g4fC  p 
mp  g2ch 
p  gSaS 
-gSaS  p 
g5a4  pd 
-g5a3  pd 


g2cb  pd 
g5a2  pd 
g5al  pd 
gSaO  pd 
mp  g3r4 
p/  g2cj 
p/  g3f3 
f  1  hov3 

gnu7 

-g3f0  p 
g4fl  pd 
-g4f2  pd 
-gdfS  pd 
-g4f4  pd 
-g4f5  pd 
~g4f6  pd 
-g4f7  pd 

gnu  8 

P  j28 
fl  k3 

m8r0 

OaOOOO 

0 

3  zz 

0 

m8rl 

080000 

0 

4  zz 

0 

m8r2 

ObOOOO 

0 

4  zz 

0 

m8r3 

080000 

0 

6  zz 

0 

m8r4 

OSOOOG 

0 

7  zz 

0 

m8r5 

OdOOOO 

0 

5  zz 

0 

m8r6 

OaSOOO 

0 

5  zz 

0 

m8r7 

080000 

0 

11  zz 

0 

m8r8 

OcOOOO 

0 

9  zz 

C 

m8r9 

08C0CC 

G 

8  zz 

m8ra 

084 COG 

0 

6  zz. 

G 

mSrb 

0c74g.C 

0 

11  zz 
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0  m9r3  0  zz 

mSrc  UafaSU  C 

0  0 

13  zz  0 

0  m88f  OLiaGGO 

mSrd  C94000  0 

0  13  zz 

G  zz  0 

0  clal  0a93fC 

mSre  OaTOOO  G 

0  16  zz 

8  zz  0 

0  mSrh  OafScO 

mSrf  080000  G 

0  15  zz 

9  zz  0 

0  mSrj  0a3c00 

mSrg  OfcOOO  0 

0  11  zz 

6  zz  0 

0  mSrk  080000 

nrlb  0  zz  0 

0  '  1C  zz 

0  0 

0  mSrn  Ce6000 

m8rt  0fedf9  C 

0  7  zz 

20  Z7  0 

0  m8rp  0b7cC0 

m8ru  088d0f  5 

cCOOOO  10  ZZ 

22  zz  0 

0  mSrq  083c00 

m8rv  0b3580  0 

Q  1C  zz 

13  zz  0 

0  mSrs  OfOSSO 

mSrw  OfOlfO  0 

0  IS  zz 

17  zz  0 

C  in9C2  GaOCOO 

m8ry  091 fSe  0 

C  4  zz 

19  zz  0 

0  iti9r4  0  zz 

m8rz  Obi 800  0 

0  0 

9  zz  0 

0  m9r5  0  zz 

m9r0  0bc200  C 

0  0 

15  zz  0 

0  ni9r6  0  zz 

m9rl  080000  0 

0  0 

5  zz  0 

0  m9r7  0  zz 

m991  OcOOOO  0 

0  c 

2  zz  0 

0  m9r8  0  zz 

m9r2  0  zz  0 

0  0 

0  0 

0  0 
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mSr'J 


mOr  a 


mOrb 


mOrc 


mOrd 


mOre 


mOrf 


mCJrg 


mOrh 


m9r  j 


m3rk 


mOrn 


m9rq 


m9rt 


mare 


mOru 


m9rw 


0  zz 
0 
0 
C 

0  zz 
0 
0 
0 

C  zz 
0 
0 
0 

0  zz 
0 
0 
0 

0  zz 
0 
0 
0 

0  zz 
0 
c 
0 

C  zz 
0 
c 
0 

0  Z7. 
0 
C 
0 

C  zz 
0 
0 
0 

0  zz 
0 
c 
c 

C  zz 
/ 

V 

0 

0 

0  zz 
0 
0 
0 

0  zz 
0 
G 
0 

0  zz 
0 
0 
G 

0  zz 
C 
G 
C 

0  zz 
0 
V 
0 

G  zz 


G 

G 

0 

mOrx  C  zz 
G 
0 
G 
C 

g2cl  GbGOGC 
G 

5  zz 
0 

g2c2  Ge4000 
G 

7  zz 
0 

g2c3  OaaCOO 
C 

9  zz 
G 

g2c4  OcSCOO 
0 

5  zz 
C 

g2c5  CdOOCO 
C 

4  zz 
0 

g2c6  OefOOO 
0 

8  zz 

0 

g2c7  OeOOOO 
G 

4  zz 
0 

g2c8  GcCOOO 
G 

3  zz 
C 

g2c9  0S1580 
C 

13  zz 
0 

g2ca  0a328G 
0 

17  zz 
C 

g2cb  0fc400 
0 

10  zz 
0 

g2cc  OeacOO 
0 

10  zz 
0 

g2cd  CeaGCO 
0 

12  zz 
C 

g2ce  OcGeGO 
0 

1 1  zz 
0 

g2cf  Ge8800 
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0 

9  zz 
0 

g2cg  OOecOO 
0 

10  zz 
C 

g2ch  Cd-IOGO 
0 

6  zz 
0 

g2cj  080000 
0 

14  zz 
0 

g2ck  G8000G 
0 

3  zz 
0 

g3f0  0  zz 
0 
0 
0 

g3fl  0  zz 
0 
0 
G 

g3f2  0  zz 
C 
0 
c 

g3f3  C  zz 
C 
0 
0 

g3f4  0  zz 
0 
0 
C 

g3  f  5  0  zz 

0 
0 
0 

g3f6  0  zz 
0 
0 
0 

g3f7  0  zz 
0 
0 
0 

g318  0  zz 

0 
0 
0 

g3f9  0  zz 
0 
0 
0 

g4f0  0  zz 


C 

G 

G 

g4fl  0  zz 
0 
C 
0 

g4f2  0  zz 
C 
C 
G 

g4f3  0  zz 
C 
0 
0 

g4f4  G  zz 
0 
c 

G 

g4f5  0  zz 
G 
0 
C 

g4fC  0  zz 
0 
0 
0 

g4f7  0  zz 
0 
0 
0 

gSaC  0  zz 
0 
G 
C 

gSal  0  zz 
0 

0 

g5a2  G  zz 
0 
0 
0 

g5a3  0  zz 
0 
0 
0 

g>,a4  0  zz 

0 
0 
0 

gSaS  0  zz 
0 
0 
0 

mO  r  f  p 

IpI 

P  g3f0 

m9rx  p 
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hov2 


hov3 


hov4 


hov5 


p 

R3f  1 

B3f0 

P 

mp 

hov2 

-R3f  1 

pd 

f  c  1 

gnu3 

f  1 

0a3d70 
a3(l7  0a 
f f f900 
000000 

k3 

P 

K4f7 

R3f0 

IpI 

P 

p 

m9r  f 

E4f7 

IpI 

P 

p 

m9rg 

mO  r  f 

P 

mp 

hov2 

-m9  r  g 

pd 

fcl 

gnu  7 

f  1 

l<3 

pc 

37 

pc 

16 

pc 

16 

pc 

16 

pc 

41 

j  126 

P 

Pt 

11 

n  1 

k6S 

j  126 

P 

mp 

mSro 

mp 

j6 

P 

j  126 

P 

•ttol 

1  s  r 
ho  v7 

hov6 

f  1 

k69 

NOTE:  Rpsorvc  in  the  memory  200  consecutive  floating  point  number  locations  for  each 
of  the  startiiiR  addresses  listed  below.  There  will  then  be  in  the  memory  2200  floating 
point  number  reservations.  For  example,  beginning  with  address  blOO  there  will  follow 
blOl,  bl02,  etc.,  each  of  which  is  reserved  for  storage. 

blOO 

clOO 

.1300 

j400 

j900 

k300 

vxlO 

vvxl 

vlOO 

wlOO 

ylOO 


♦  ♦  ♦ 
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REPRODUCTION  OF  PROGRAM  IN  MACHINE  LANGUAGE  FORM 

The  following  list  of  commands  is  a  Narec  “machine  language”  translation  of  the 
symbolic  program  appearing  in  Appendix  B.  It  is  to  be  noted  that  list  contains  both  fixed- 
point  and  floating-point  commands.  Two  columns  of  numbers  appear  in  the  list.  The 
left  hand  column  is  extraneous  to  the  actual  computation  schedule.  It  is  used  only  for 
visual  identification  of  the  lines  in  the  right  hand  column  which  constitutes  the  true  read- 
in.  The  left-hand  column  is  not  read  into  the  computer.  However  the  lines  in  the  program 
must  be  read  in  consecutively  beginning  at  address  0400. 

Each  line  of  command  consists  of  twelve  hexadecimal  digits.  They  are  interpreted 
by  the  computer  logic  from  left  to  right  as  follows; 

digits  1  through  4  represent  one  address  in  memory 

digits  5  and  6  represent  a  fixed-point  or  floating  point  order 

digits  7  through  10  represent  one  address  in  memory 

digits  11  and  12  represent  a  fixed-point  or  gloating-point  order. 

Note  that  the  program  starts  in  fixed-point.  To  enter  floating-point  as  a  particular 
address  one  uses  the  standard  form 

(address)  50  0200  11 

To  exit  floating  point  into  fixed  point  one  uses  the  symbols  E4  or  E5. 

The  symbol  H  indicates  the  division  of  a  Narec  word  between  the  left  order  and  the 
right  order. 

The  symbol  *  means  “skip  this  order.” 
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TYPICAL  OUTPUT  SHEET 


A  typical  output  sheet  is  reproduced  below  with  added  description  of  the  meaning  of  the 
numbers.  A  number  without  a  decimal  point  is  interpreted  as  fraction  followed  by  an 
exponent.  For  example 


is 


+6S511109498  +02 

+65.511109  498 


29 
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nnnn 


0000 

■  +0.000000000  -2 

I  +6. 

h  +6.60O 


0000 

0001 

+65511109498 

+02 

0000 

+10000000827 

-07 

0002 

+10936175364 

-01 

0003 

+96013725118 

-01 

0004 

+38634237703 

+00 

0005 

+83439671453 

+00 

0006 

-39362130016 

+00 

0007 

+73217232273 

-01 

0008 

-78111671012 

-02 

0009 

+55342208234 

-03 

0010 

-2624o348o44 

-o4 

0011 

+10937705008 

-05 

0012 

-33366009971 

-07 

0013 

+82430513625 

-09 

0014 

-16852080723 

-10 

0015 

+29012288326 

-12 

0016 

-42668929896 

-l4 

0017 

+54260160987 

-16 

0018 

-60277682314 

-18 

0019 

+59019532925 

-20 

0020 

-51329455278 

-22 

0021 

+39924268608 

-24 

0022 

-27941066789 

-26 

0023 

+17690620699 

-20 

0024 

-10102634883 

-30 

0025 

+53519977323 

-33 

0026 

-25790419362 

-35 

0027 

+11436205148 

-37 

0020 

-468232U649 

-4o 

0029 

+17755655572 

-42 

0030 

-62540105973 

-45 

0031 

+20515535010 

-47 

0032 

-62832443777 

-50 

0033 

+10007915429 

-52 

0034 

-484oi3o6663 

-55 

0035 

+12224704706 

-57 

0036 

-29069549227 

-60 

0037 

+65196116625 

-63 

0038 

-13814023500 

-65 

0039 

+27696270773 

-60 

oo4o 

-52622999510 

-71 

004l 

+94805173204 

-74 

0042 

-16258237917 

-76 

0043 

0044 

+26506679241 

-41169066792 

Oo45 

+60904915704 

-05 

0046 

-86255720730 

-80 

0047 

+11660682339 

-90 

oo40 

-15082269790 

-93 

0049 

+10602422465 

-96 

0050 

-22103121071 

-99 

Ei  genva 1 ue 

Precision  of  Eigenvalue 


Prolate  Spheroidal  Wave 
Expansion  Constants 


APPENDIX  E 
NOTES 


The  computation  results  shown  in  Appendix  D  requires  116  seconds  for  completion. 
This  time  depends  largely  on  the  precision  of  the  eigenvalue  deemed  acceptable.  The 
precision  is  controlled  in  turn  by  the  number  of  times  the  corrections  is  applied  to 
the  eigenvalue  a^,  •  In  this  program  the  number  of  corrections  is  set  at  four.  If  greater 
precision  is  demanded,  one  enters  the  program  on  the  right  side  of  location  069C  and 
replaces  the  existing  address  (=058 F)  containing  the  existing  number  of  corrections 
(minue  one)  with  another  address  containing  a  larger  number.  For  example  if  it  is  desired 
to  apply  9  corrections  to  instead  of  4  corrections  one  scans  through  the  program,  to 
find  the  fixed  number  “8."  This  is  seen  to  be  at  address  0793.  Hence  location  069C  is 
changed  to  read 


060F  F8  0793  DC 

The  computation  of  A^,!  is  highly  sensitive  to  the  starting  value  inserted  in  address 
05B9.  When  h  is  large  (i.  e.  >10)  we  use  a  formula  by  Meixner  (5)  to  help  fix  a  good 
approximation  for  A^, .  The  Meixner  approximation  has  eight  terms  (prolate)  or  seven 
terms  (oblate).  In  the  program  the  first  term  in  either  case  is  multiplied  by  a  small 
number  (io"^/2)  fixed  in  location  086B  and  086C.  The  last  term  is  then  compared  with  the 
resultant  multiplication.  If  the  last  term  is  l/lOO  or  less  of  the  first  term  then  the 
Meixner  approximation  is  used.  Otherwise  the  standard  approximation  (i.e.  approximation 
h  -  0)  is  used. 
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0400 

063810000000 

0401 

040150020011 

0402 

0585e60583ea 

0403 

0593e80585e6 

0404 

058be90595e8 

0405 

0585e6058deb 

0406 

0597e8058bea 

0407 

0599e80591e9 

0408 

059be80583e6 

0409 

0583eb058deb 

040a 

059de8058deb 

040h 

059fe80583e6 

040c 

058be905ale8 

040d 

0583e6058deb 

040f 

05a3e8058bea 

040f 

05a5e8058de9 

0410 

05a7e80591ea 

0411 

05a9e80591e9 

0412 

058de905abe8 

0413 

0587e60587eb 

0414 

05ade8058dec 

0415 

05af e805ade6 

0416 

05adeb05ble8 

0417 

0585e60595eb 

0418 

05b3e80597e6 

0419 

0595eb059dea 

041a 

058bea05b5e8 

041b 

0599e6059beb 

041c 

05b7e805b5e6 

041d 

05b7ec05adeb 

041e 

05b3e906d4e9- 

041f 

0420 

0421 

05b9e806d6e2 

0593e60^bbe8 

05bbe6058bdc 

Go  to  left  order  of  0638  |1  * 
Enter  floatingpoint 


Compute  “Starting*  = 
r 


2^(-C-ti)  -  2m^  - 
(2^1-1)  (2\+  3) 


Store  in  05B9  ||  Go  to  left  order  06D6 
Is  -t  -m  an  odd  number  or  an  even  number? 


l.OOOO _ 050 
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0422 

0426u20428c2 

0423 

0424e20000e0 

0424 

05bbe6058dea 

0425 

05bbe80421e2 

0426 

05bde605bf e8 

0427 

05cle80445e2 

0428 

058be605b£e8 

0429 

05cle80445e2 

042a 

0583e605cle9-^ 

042b 

05c3e8058be9 

042c 

05c3eb05c5e8 

042d 

059fe6058bea 

042e 

05c7e805cle6 

042f 

058deb05c9e8 

0430 

05a5e605c9e9 

0431 

05cbe80591e9  * 

0432 

05cbeb05cdo8 

0433 

05c7e605cdec 

0434 

05cfe8058be6 

0435 

05cf ea05af eb 

0436 

05c5e905dle8 

0437 

0000e30000e0 

0438 

05cle6058bea 

0439 

05d3e805a3e6 

043a 

05cle905d5e8 

043b 

058bea05d7e8 

043  c 

05cle905d9e8 

043d 

058de905dbe8 

043e 

0591ea05dee8 

043f 

05cle605d3eb 

0440 

05d5eb05d7eb 

0441 

05d9ec05d9ec 

0442 

05deec05dbec 

0443 

05bleb05e0e8 

0444 

OOOOeSOOOOeO^ 

( -t  -  is  even)  bring  up  zero  1|  Store  zeroin  05BF 
Store  zero  in  05C1  ||  Go  to  0445 


(-t  -  m  is  odd)  Bring  up  *1"  jj  Store  "I*  in  05BF 
Store  ‘1*  in  05C1  H  Go  to  0445 


Compute  T®  =  (m  +  n)  (m  +  n  +  1) 


"T 


_ 4m^-  1 _ 

(2m  +  2n  -  1)  (2m  +  2n  +  3) 
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.0445 
0446 
0447 
0448 
0449 
044a 
044b 
044c 
044d 
044  e 
044f 
0450 
0451 
04  52 
0453 
0454 
0455 
04  56 
0457 
0458 
0459 
045  a 
045b 
045c 
045cl 
045e 
045f 
0460 
0461 
0462 
04  63 
0464 
0465 
0466 
0467 

0468 

0169 


046  a 
046b 
046c 
040d 
046e 
04  6  £ 
0470 
0471 
0472 
0473 
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058bc61900e8 

0593e60591e9 

05e2e805bfe6 

05e4e805cle8 

05e4e6058de9 

05e4e805e2dc 

044de20456e2 

0456e20000e0 

044ee40000e0 

045050044110 

044420044422 

045050020010 

0438e205e0e6 

0890c80000e0 

0002d0045200 

05cle6058de9 

05cle80449e2 

05bf e605cle8 

05e6e80000e0  ,, 

05e6e6058de9 

05e6e805e2dc 

045ce20465e2 

0465e20000e0 

045de40000e0 

045f50045el0  i 

043720043722  | 

045£50020010  , 

042ae205  dle6 

0a20e80000e0 

0002d0046100 

05cle6058de9 

05cle80458e2 

05bf e605cle8 

05e8e80000e0 

05ebe6058de9 


Compute  [Nn  +  j]^ 


+  A. 


T 


m 

n 


05o8o805ij2dc 

046be20474o2 

0474e20000e0 

0890c70000o0 

lOOOecOOOOcO 

0a20ea05b9t'9 

1902c805eae8 

0002d0046b00 

0002d0040c00 

0002d0046d0n 


Fi ni shed? 

No:  Go  to  046B  II  Yes:  Go  to  0474 
Yes:  Go  to  0474  |1  ♦ 


II  Store  in  05EA 


0002d0046e00 
087FE200O0E0  ^ 
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0475 

059 je6058de9-^ 
05ece805f2eS 

0476 

0477 

0478 

0479 

058dea049feH 
05bfe605cee8 
05eee605  8de9 
05eee805  r2dc 

047  a 

047cc20485ol. 

047b 

0485e20<''f'f'.'() 

047c 

047d 

047e 

047f 

0480 

0481 

0482 

0483 

0484 

0892e60000e0 

1902ec0000c0 

1902ec0000c0 

ObbOeSOOOOcO 

0002d0047c00 

0002d0047d00 

0002d0047f00 

0002d0047e00 

0478e20000e0— 

0485 

047fd8048700-^ 

0486 

058be60000f0 

0487 

0000u80000e0 

0488 

0bb0u60000c0 

0489 

058bdc048bc2 

048a 

04alc2048bf2 

048b 

0488d8049100 

048c 

0488d8048e00 

048d 

0002d0049100 

048e 

0000e60d40e8 

048f 

05bfe605f0e8 

0490 

0d40e60000e0 

0491 

OOOOebOOOOeO 

0492 

0d40e80000c0 

0493 

05f0e6058do9 

0494 

0495 

0496 

0497 

0498 

0499 

049a 

049b 

049c 

04  9d 

04  9e 
049f 
04a0 

05f0e8049fdc 

0497e20499e2 

0499e20000e0 

0002d0049100 

0490e20000e0 

0002d0048800 

0002d0049200 

0002dl048e00 

0002d0049000 

049fe6058dea 

049fe80488e2 

000000000000 

000000000000^ 

Compute 


for  X  =  1,  2,  .  •  • 


II  Finished? 

No:  Go  to  047C  ||  Yes:  Go  to  0485 

Yes:  Go  to  0485 


Finished? 

Yes:  Go  to  04A1 


No:  Go  to  048B  ^ 
No:  Go  to  048B  J 


— U  Outer  Loop 


Compute 


> 

1 

3 

1 

etc. 


II  Fini,shed? 

No:  Go  to  0497 
Yes:  Go  to  0499 


II  Yes: 
11* 


1  Inner 

Go  to  0499  /Loop 
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O-iri2d904a200-*- 


05Sbo60000e8 


0dl0o6058bdc 


0du6e204abe2 

04a6c20000e0 


Compu  t  c 


, ,  n-  ^  n-~  , 


n5f 4e60000e0 


0d40e905f4e8 


Store  computation  in  05F4* 


0002d004a300 


0002d004a700 


04a3c20000e0 

05f4e6058be9 


05f4e80000e0- 


0593e605f6e8- 


05f6e6058de9 


05f6e805f8e8 


05f ae805cle8 
04b2o40000e0 
04b45004b310 
044420044422 
04b4G0020010 
0438e204b6e4 


Storemax.  (n  +  2x)  in  05F8 


04b85004b710 

043720043722 


-Compute  for  determining  (n  2x)  max. 


04b850020010 


042ae205dle6 


05b9t;a05fce8 


05feub0000f9 


0600e805e0e7 


05fcec0000f9 


OCOOdc04cOe2 


04c0e204aee2* 


Finished? 


lives:  Go  to  04CO 


Yes:  Go  to  04CO  H  No:  Go  to  04AE 


NAVAL  RESEARCH  LABORATORY 


37 


04c0 

04cl 

04c2 

04c3 

04c  4 

04  c  5 

04c6 

05bfc605cle8- — 

0593c6058doa 

0602c80000o0 

05f8fc'6058dea 

05f8e80000e0 

0602cG05f8oa 

04dbf40000c0 

Iff  ini  shed  go  to  04DB  ||  * 

04c7 

05f ao605cle8 

04c8 

04c9c40000c0 

04c9 

04cb5004cal0 

04ca 

044420044422 

04cb 

04cb50020010 

04cc 

0438c204cde4 

04cd 

04cf5004cel0  * 

1 -  Compute  - jj - ±L  for  x=  1,2,... 

+m  +  2x 

04ce 

043720043722 

04cf 

04cf50020010 

04d0 

042ae205dle6 

04dl 

05b9ea0604c9 

04d2 

0606e805e0e7 

04d3 

0606ec0604e8 

11  Store  +  in  0604 

04d4 

0604e60604eb 

04d5 

OSeOecOOOOeO 

04d6 

OedOeSOOOOcO 

Store  +  in  OEDO  ||  *. 

04d7 

05fae6058dea 

04d8 

05fae80000e0 

04d9 

0002d004d600 

04da 

04c3e20000e0, _ 
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04db 

04d6d904dc00^ 

0<ldc 

05SbeG0000c8 

04dd 

05Sbc60608c8 

04de 

0ed0cG058bdc 

Oldr 

04clc204eXe2 

04e0 

0  h  lc20000e0 

04ol 

OcdOoGOcSbdc 

01o2 

0ic1c204cae2 

04  0  3 

04clc20000c0 

Ole  4 

OGOScGOOOOeO 

04  c  5 

OcdOcbOOOOeO 

04o6 

10C0c80C08e8  < 

04c7 

0002d004ol00 

+  1  .  . 

+  4 

04c8 

0002d004c500 

04cn 

04elo20000c0 

04c‘a 

0002d004dc00 

04  cb 

04dodo04el00 

04cc 

04dcdS04c500 

C'lcd 

C002d004cG00 

04cc 

04ddc20000c0 

04c  I 

04cGdC04X000 

04f0 

058bcGOOOOe8 

04X1 

1060cG058bdc 

1  IlFinished? 

04X2 

04X4e204X9e2 

No:  Go  to  04F4  ||  Yes:  Go  to  04F'> 

04X3 

04f4c20000c0 

No :  Go  t o  04F4 

04X4 

OGOacGOOOOcO 

11  Store  y  + 

.1?’  1 
t.  -  m  +  4-' 

04X5 

lOGOcOOGOacf; 

et  c . in  060A* 

04X6 

0002d004X100 

04X7 

0002d004X500 

04X8 

0  ]Xlc2U(;GCe0  — 
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04  f<) 

05f 'icGUtiOucU 

BriHRup  1  +  i  etc  11  Add  2  (N|_„  ,  ^  +  2 

04  lu 

OG0cufc06O-i  L‘0 

Store  in  060C  ||  Bringup  [N|'_^^2^1 

04  fb 

05cuc'a060ccc 

Suljtract  (Nl^_„^2l2  II  Divide  by  060C  to  form 

04fc 

060ccS05bO('0 

Store  bA^^  in  060E  ||  Add  A^^ 

04fd 

05b9c80G0c(jG 

Store  -f  in  05B9  ||  Bring  up  bA^^ 

04i'c 

OOOOl'JUGlUcu 

Form  abs  value  ||  Subtract  tolerance  in  0610 

04ff 

0GlbI40Gc5t2 

Does  A|^^meet  precision?  No;  Go  to  0693  |1  Yes;  Go  to  06C5 

0500 

00Cbl'00S73L-2 

Output:  Print  to  11  significant  digits  H  Go  to  0873 

0501 

050111050210 

Go  to  right  side  of  this  line  ||  Go  to  left  order  of  0509 

0502 

0G32i.-l>05(jlcb 

Bring  up  zero  from  0632  H  Store  zero  in  05E4 

0503 

05i.-Gt.805(.-bcF 

i  n  05E6  ||  ~  i  n  05E8 

030-1 

C5(jet805  r0(j8 

-“in  OSEE  11  ~  in  05F0 

0505 

05f 2e80G0afb 

~  in  05F2  H  ~  in  060A 

0506 

05X4tjb05fGub 

~in  05F4  11  ~  in  05F6 

0507 

OoXbuSOOfucb 

~  in  05F8  1|  ~  in  05FA 

05  OS 

0604(.'SOG47(j4 

~  in  0604  11  ~Enter  fixed  point  at  0647 

0509 

050950G20C11 

Enter  floatingpoint 

05Ca 

05bboG1900(_b*^ 

Bring  up  “1"  H  Store  in  1900 

050b 

OGGccGOOCGcO 

Dummy  command 

050c 

050bUS04GaOO 

Change  exit  from  046A.to  066E 

OoOd 

0593>jG0591i,'b 

Bring  up  n-m  |1  Add  “4” 

050'. 

050f 

05c'2ob0G32cG 

05c4cS05(;Gt;8 

Store  n-m■^4m05E2  !1  Bring  up  zero 

Store  zero  in  0524  11  Store  z,  ero  in  05  E6 

0510 

05cSo80000u0 

Store  zero  in  05E8  ||  * 

0511 

063408045200 

0512 

0633cl804G100 

0513 

0634080  IGbOO 

j  f'ompiitc'nMfNy  ] 

*  ’  *■  1  ”  m  ' 

0514 

063608046c00 

0515 

0635O8O4GO0O 

1 

i 

0516 

003708046u00 

i 

0517 

0447e30000e0-^ 
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0ol8 
051 P 
051  a 
051b 
051c 
rjia 
051c 
0511 
0520 
0521 
0522 
0523 
0521 
0525 
0526 
0527 
0528 
0529 
052  a 
052b 
052c 
0520 
052e 


069209051900 
058bcCOOOOc8 
05blo00612o3 
0012c60580c0 
0612eSOOOCcO 
1902c60000X9 
053b0c0520o2 
052fc2052Uo2 
0612c6058dcb 
05a3c9G5ybo9 
0611c8058dea 
OGl  Icb0616o8 
0612cG03a3c9 
0618oS058bca 
0618eb061ao8 
0616e6061aoc 


Transfer  19xx  to  right  of  0519 

Bringup  “1"  ||  Store  “1"  in  1902 
Bringup  zero  (even)  or  1  (odd)  1|  Store  in  0612 
Bring  up  zero  or  “1"  ||  Add  2 
StoreO(orl)  +  2x  in  0612  (  =  counter)  ||  * 

Bring  up  1902  +  2x  ||  Form  abs  value 
Finished?  |1  No;  Go  to  0520 
Yes;  Go  to  052F  |1  No;  Go  to  0520 
BringupO(orl)  +  2x  H  Mult,  by  2 
Add  2m  ||  add  1 
Store  in  0614  ||  subtract  2 
Multiply  by  0614  |!  Store  in  0616 

BringupOfor  1)  12  x  |1  add  2m 
Store  in  0618  ||  Subtract  1 

Multiply  by  0618  |1  Store  in  061A 
Bring  up  0616  ||  Divide  by  061A 


05adoc0000o0 

1902cb0000e0 

1110o80000c0 

0002d0051d00 

0002d0032900 

0002d0052a00 

051be20000e0 


Comput  e  al 1 
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0521 

0530 

0531 

0532 

0533 

0534 

0535 

0536 

0537 

0538 

0539 

053a 

053b 

053c 


052ad9053000' 

058be60000e8 

061ce6058be9 

061ce80000e0 

Ilf0e60000f9 

058bdc0536c2 

Q53de20536e2 


Compute  '  “T" 


Finished?  II  No:  Go  to  0536 
YesT”Go  to  053D  ]]  No:"  Goto  0536 


061ce60000e0 

llfOebOOOOeO 

1380e8061ce8 

0002d0053300 

0002d0053700 


0002d0053800 

0533e2000r'c0-.- 


e  t  c  . 


053d 

053e 

053f 

0540 

0541 

0542 

0543 

0544 

0545 

054G 

0547 

0548 

0549 

054a 

054b 

054c 

054d 

054e 

054f 

0550 

0551 

0552 

0553 


OSOadOOuGcOO-^ 

058bu600C0c8 

05bXuC061eo8 

OGleuOOubdc'O 

061ec‘8061cc6 

05a3u90620cS 

0622e80000e0 


Compute  ^  ^ 

n  !  d  _ 


1380e00000f9 

058bdc0547e2 

0555c20547o2 

061ce60000f9 

0620dc054af2 

054ee2054c(j2 

OC20<i6058bea 

0624e80622ob 

0622eS0624c6 

0620e80547e2 

0622e6061eec 

1380eb0000a0 

ISlOeSOOOCfO 

0002d0054400 

0002d0054X00 

0002d0055000 


Finished?  Go  to  0555 


0554 


0540e20000e0 


42 


NAVAL  RESEARCH  LABORATORY 


05G5 

0550d9055600 

0556 

05RbcG0000c8 

0557 

1510c60000f9 

0558 

05P.bdc055ac2 

0559 

065ce2055ae2 

055a 

0C2GcG0000o0 

055b 

1510o90626c8 

055c: 

C002d0055700 

055d 

0002d0055b00 

055e 

0557e20593c6. 

055f 

0632UcU665c'3 

0560 

056oc3062ac8 

0561 

05Sbca062cc8 

0562 

0632dc0566o2 

0563 

0566o2062ac6 

0564 

062ccb062aoS 

0565 

062co60561c2 

0566 

0585c60583e9 

0567 

0632dc0665c3 

0568 

0570c20000o0 

0569 

062cc8008boa 

056a 

062bc80632dc 

056b 

0571c20571c2 

056c 

062co60628eb 

056d 

062cc80628c6 

056c 

0569o305Sbc6 

0561 

062ac80566c2 

0570 

05Sbc6062cco 

0571 

062cc6yG2ucc 

0572 

0626cc0630cS 

0573 

IGaOcSOOOOcO. 

Compute  y  +  2m)! 
Z__i  n ; 


Finished?  {|  No;  Go  to  OSSA 
Yes:  Go  to  06SC  II  No:  Go  to  OSSA 


J  Return  to  05S7  ||  Bring  up  ‘-t-m" 

-  Compare  with  zero  ||  (-t-m)  <0.  Error:  Go  to  right  of  066 S 
(-C-m)  =  0.  Go  to  right  of  0S6E  ||  ("t-m  )  >  0.  Store  (-i'm)  in  062A 
Subtract  “I*  ||  Store  in  062E 

Compare  with  zero  ||  (-t'm-1)  <  0.  Go  to  left  side  of  0S66 
(-t-m-1)  =  0.  Gotoleftsideof0S66||(^-m-l)>0.  Bringup  062A 
Mu  1 1 i ply  by  062E  |1  Store  in  062A 
Bring  up  06  2E  ||  GotoleftsideofOS61 
Bring  up  ||  arid  A^^ 

Compare  wi  th  zero  ||  ( -{-m)  <0.  Error:  Go  to  right  06  6S 

- Compute  dfi  or  dj  ('t+m)  =  0.  Go  to  left  of  0S70  ||  * 

(^+m)  >  0.  Store  (-t+m)  in  062C  ||  Subtract  “1* 

Store  (^+m-l)  in  0628  1|  Compare  with  zero 

<  0.  Go  to  lef  t  of  0S71  1|  (-t+m-l)  =  0.  Go  to  left  of  0S71 
(•t+m-1)  >  0.  Bring  up  062C  H  Mu  1 1  i  ply  by  0628 
Store  resul  t  in  062C  1|  Bring  up  0628 
Go  to  right  side  of  0S69  11  Bring  up  1 
Store  in  062A  l|  Go  to  loft  of  0S66 
Bring  up  "1"  ||  Store  in  062C 
Bring  up  062C  II  Di  vi  de  l)y  062A 
Divide  by  0626  II  Store  d^^  ^  in  0630 

-  Store  d^  j  in  16A0  I1  * 
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0574  138OcUOOO0f9 
0575  058bac0577u2 
0576  057dc'20j77c2 
0577  1380fG0C30cb 
0578  16a2e80000e0 
0579  000260057400 
057a  000260057700 
057b  000260057800 
057c  0574o20000c0 
0576  00066oU64ac5 
057e  000000000000 
057f  16a0c60000c0 
0580  OOOOcOOOOOcO 
0581  000260057X00 
0582  057<jo30000e0 
0583  080000000000 
0584  000100000000 
0585  080000000000 
0586  000200000000 
0587  080000000000 
0588  000100000000 
0589  OcOOOOOOOOOO 
058a  000400000000 
058b  080000000000 
058c  000100000000 
0586  080000000000 
058e  000200000000 
058f  OcOOOOOOOOOO 
0590  000200000000 
0591  080000000000 


0592 


000300000000 
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05b3 

000000000000 

0ob4 

000000000000 

051)5 

000000000000 

05bG 

000000000000 

05b7 

000000000000 

05b8 

000000000000 

05b9 

000000000000 

05ba 

000000000000 

03bb 

000000000000 

05bc 

000000000000 

05bU 

000000000000 

05bc’ 

000000000000 

05bf 

000000000000 

05c0 

000000000000 

05cl 

000000000000 

05c2 

000000000000 

05cu 

000000000000 

05c4 

000000000000 

05c5 

000000000000 

05c6 

000000000000 

05c7 

000000000000 

05c8 

000000000000 

05c9 

000000000000 

05ca 

000000000000 

05cb 

000000000000 

05cc 

000000000000 

05cd 

000000000000 

05ce 

000000000000 

05cf 

000000000000 

05d0 

000000000000 

05dl 

000000000000 

05d2 

000000000000 
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05d3 

000000000000 

05d4 

000000000000 

05d5 

000000000000 

05d6 

000000000000 

05d7 

000000000000 

05d8 

000000000000 

05d9 

000000000000 

05da 

000000000000 

05db 

000000000000 

05dc 

000000000000 

OSdd 

OOOOOOOSdelO 

05de 

000000000000 

03d  f 

000000000000 

05e0 

000000000000 

05el 

000000000000 

05e2 

000000000000 

05e3 

000000000000 

05e4 

000000000000 

03e5 

000000000000 

05e6 

000000000000 

05e7 

000000000000 

05e8 

000000000000 

05e9 

000000000000 

05ea 

000000000000 

05eb 

000000000000 

OSec 

000000000000 

03ed 

000000000000 

05ee 

000000000000 

05ef 

000000000000 

05f0 

000000000000 

Ojfl 

gooooooooooo 

05f  2 

000000000000 
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05f3 

000000000000 

05f4 

000000000000 

05f5 

000000000000 

05f6 

000000000000 

05f7 

000000000000 

05fS 

000000000000 

05f9 

000000000000 

05fa 

000000000000 

05fb 

000000000000 

05f  c 

000000000000 

05fd 

000000000000 

05f  e 

0a3d70a3d7a0 

03ff 

fffaOOOOOOOO 

0600 

000000000000 

0601 

000000000000 

0602 

000000000000 

0603 

000000000000 

0604 

000000000000 

0605 

000000000000 

0606 

000000000000 

0607 

000000000000 

0608 

000000000000 

0609 

000000000000 

060a 

000000000000 

060b 

000000000000 

060c 

000000000000 

060d 

000000000000 

060e 

000000000000 

060f 

000000000000 

0610 

OafecOOOOOOO 

0611 

ffdcOOOOOOOO 

0612 

000000000000 
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0613  000000000000 

0614  000000000000 

0615  000000000000 

0616  000000000000 

0617  000000000000 

0618  000000000000 

0619  000000000000 

061a  000000000000 

061b  000000000000 

061c  000000000000 

061d  000000000000 

061e  000000000000 

061f  000000000000 

0620  000000000000 

0621  000000000000 

0622  000000000000 

0623  000000000000 

0624  000000000000 

0625  000000000000 

0626  000000000000 

0627  000000000000 

0628  000000000000 

0629  000000000000 

062a  000000000000 

062b  000000000000 

062c  000000000000 

062d  000000000000 

062e  000000000000 

062f  000000000000 

0630  000000000000 


0631 


000000000000 
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0632 

000000000000 

0633 

000000000000 

0634 

089000000000 

0635 

0a2000000000 

0636 

190000000000 

0637 

190200000000 

0638 

063d32063910 

0639 

063d81058308 

063a 

064581057d01 

063  b 

040083240030 

063  c 

040110063dl0 

063d 

000000000000 

063e 

000000000000 

063f 

000000000000 

0640 

000000000000 

0641 

oooooooooooo 

0642 

000000000000 

0643 

oooooooooooo 

0644 

oooooooooooo 

0645 

oooooooooooo 

0646 

oooooooooooo 

0647 

2400810400fd 

0648 

064850020011 

0649 

05bfe605cle8 

064a 

0445e2000082 

064b 

240083040030 

064c 

063810000000 

064d 

002513002513 

064e 

000713000413 

064f 

0004130S83e6 

Read  m,  -t,  h,  etc.  into  063D  etc.  ||  Go  to  0639 
Transfer  8  data  lines  from  06 3D  etc  into  0583  etc. 
Transfer  1  data  line  from  0645  into  057D 
Transfer  768  commands  from  0400  etc.  to  2400  etc. 
Go  to  left  order  0401  ||  * 


Transfer  fresh  program  (256  lines)  from  2400  to  400 
Enter  floatingpoint 

Bring  up  “odd-even"  (=  “one"  or  “zero")  1|  Store  in  OSCl 
Go  to  0445  II  COMPUTATION  FINISHED:  STOP 

Output:  Carriage  return  (C.R.)  ||  C.R. 

Print  character  "m*  ||  Space  (on  typewriter) 

Space  11  Bring  up  “m" 
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0650 

0300fl0025f3 

Print  value  of  m  ||  C  R- 

0651 

0006f30004f3 

Print  character  “-t"  ||  Space 

0652 

0004f30535e6 

Space  II  Bring  up 

0653 

0300fl0025f3 

Print  value  of  “-t*  ||  C.  R. 

0654 

0008f30004f3 

Pr  i  nt  cha  rac  t  e  r  “  h  "  ||  Space 

0655 

0004£30587e6 

Space  II  Bring  up  “h* 

0656 

0303fl0025f3 

Print  value  of  “h"  ||  C.  R. 

0657 

0025f30025f3 

C.R.  II  C.R. 

0658 

0010f30010f3 

Print  character  “0*  ||  “0" 

0659 

0010f30010f3 

“0"  II  “0* 

065  a 

0029f305b9e6 

Tab  (on  typewr  i  t  er)  ||  Bring  up 

065b 

0500e20000e0 

Go  to  0500  II  ♦ 

065c 

05a3e60632dc 

Bringup2m  ||  Compa  re  2m  with  zero 

065d 

0665e30666e3 

2m  is  <  zero.  This  is  an  error.  Go  to  right  0665  ||  2m  =  0 ,  Go  to  0666 

065e 

066ce2058bea 

2m  >0.  Gotoleftof  066C  ||  Subtract  “1’ 

065f 

066ae80632dc 

Store  066A  ||  Compa  re  w  i  th  z  o  ro 

0660 

0663e30663e3 

Result  is  less  than  zero.  Go  to  right  of  0663  ||  Result  equals 
zero,  eo  to  right  of  0663 

0661 

06C8e6066aeb 

Result  >0.  Bring  up  0668  ||  Mu  1 1  iply  by  066A 

0662 

0668e8066ae6 

Store  in  0668  ||  Bring  up  066A 

0663 

065ee30626e6 

Go  to  right  of  065E  ||  Bring  up  0626  (^working  cell) 

0664 

0668eg0626e8 

Add  number  in  halo  1|  S  t  or  e  r  esu  1  t  i  n  0626 

0665 

055ee30666e4 

Go  to  right  of  05 5E  ||  Enter  fixed  point  at  0666 

0666 

0bad82058be6 

Erroi  stop  ||  Bring  up  “1" 

0667 

0668e80663e3 

Store  "1"  in  0668  ||  Go  to  right  side  0663 

0668 

000000000000 

0669 

000000000000 

066a 

000000000000 

066b 

000000000000 

066c 

05a3e605bfe9 

Bring  up  2m  |1  Add  “0"or  “1"  f rom  05BF 

066d 

0668e8065ee3 

Store  in  0668  |I  Go  to  right  side  of  065E 
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066  e 

05b£e605cle8 

Bring  up  “0*  or  “1”  1|  Store  “0”  or  “1”  in  05C1 

066f 

0589e603f8e8  ' 

Bring  up  upper  limit  of  constants  jj  Store  in  05F8 

0670 

05£aey0632e6 

Store  in  05FA  ||  Bring  up  zero 

0671 

0604e8058be6 

Store  “0"  in  0604  ||  Bring  up  “1" 

0672 

lE02E80593Er) 

Store  “1"  inlE02  ||  Bringup-t-m 

0673 

058DEA0602E8 

Subtract  “2*  ||  Store  -t-m-2  in  0602 

0674 

05f8e6058dea 

Bring  up  05F8  II  Subtract  2 

0675 

05f8o80000e0 

Store  in  limit  cell  (-  05F8)  il  ♦ 

Q676 

0602e605f  8tja 

Bringup  0602  II  Subtract  (limit 

-  2x) 

0677 

06a8£40000e0 

If  0602  <  (limit  -2x)  Go  to  0678  (left  side):  If  0602  > 

(1  imit  -2x)  Go  to  06AB||  * 

0678 

05£ae605cle8 

0670 

067ac40000c‘0 

067  a 

067c50067bl0 

^  Final  calc,  of  +  ^ 

}.  Calc. 

067b 

044420044422 

to  Nj  (k  =  max.  subscript)  ' 

067  c 

067C50020010 

0676 

0438e2067ee4 

✓ 

067e 

068050067£10 

1 

067  f 

043720043722 

1 

0680 

068050020010 

1 

>-  Calc.  T 

0681 

042ae205(lle6 

J 
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0682 
0682 
0684 
0685 
0686 
0687 
0688 
0689 
068a 
068b 
U68c 
U68d 
068e 
068  f 
0690 
0691 
0692 
0693 
0694 
0695 
0696 
0697 
0698 
0699 
069a 


058be60000e8 

0518e20000e0-«i  li  * 
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0696 

069£e6058be9 

Br  ing  up  counter  base  (=  zero)  ||  Add  one 

069  c 

069fe8058fdc 

Store  count  in  069F  ||  Is  this  number  greater  than  3? 

069d 

0502e20502e2 

Number  less  than  3.  Go  to  0502  ||  Number  equal  to  3:  Go  to  0502 

069e 

06ale20000e0 

Number  greater  than  “3"  .  Go  to  06A1  1|  ♦ 

069f 

OOOOOOOOOOOO 

06a0 

000000000000 

06al 

05f ee606a4ec 

Br  ing  up  10“^/x  ||  Divide  by  “2’  (=  06A4) 

06a2 

05£ee80632e6 

Store  10“^/x(x  =  l,  2,  4,  etc)  in  05FE  ||  Bring  up  zero 

06a3 

069fe8087ae2 

Store  in  069F  ||  Go  to  087A 

06a4 

080000000000  ] 

06a5 

000200000000  J 

2 

06a6 

OOOOOOOOOOOO 

06a7 

OOOOOOOOOOOO 

06a8 

06b5e60000£9-«- 

06a9 

06b7e80000e0 

06aa 

06b0d8067700 

06ab 

06bldH067200 

Increase  “n"  in  (d^,),^^  by  10.  Then  recompute 
all  N^  (k  =  limit  -2x) 

06ac 

06b2d8068500 

06ad 

0589e606a4e9 

06ae 

06b3e805f8e8 

06af 

0670e20000e0-<- 

06b  0 

06b900000000 

06bl 

le0200000000 

06b  2 

leOOOOOOOOOO 

06b3 

OOOOOOOOOOOO 

06b4 

OOOOOOOOOOOO 

06b5 

OOOOOOOOOOOO 

06b6 

OOOOOOOOOOOO 

06b7 

OOOOOOOOOOOO 

06b8 

OOOOOOOOOOOO 
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06bg 

06b7e706b5fd 

06ba 

0000f906c3dc 

06bb 

0689e20689e2 

06bc 

06b5e60000f9 

06bd 

06b7e80000e0 

06be 

06bld8067200 

06bf 

06b2d8068500 

06c0 

06b3e606a4e9 

06cl 

06b3e805f8e8 

06c  2 

0670e20000e0 

06c3 

OafecOOOOOOO 

06c4 

ffdcOOOOOOOO 

06c5 

05b9e60000f9 

06c6 

06d2e80000e0 

06c7 

06c9d904ff00 

06c  6 

069be20000e0 

06c  9 

ObcaOOOOOOOO 

06c  a 

05b9e60000f9 

06cb 

06d2ea0000f9 

06cc 

06d0ea06cef 4 

06c  d 

064de2000000 

06ce 

05b9e60000f9 

06c  f 

06d2e8069be2 

06d0 

OafecOOOOOOO 

06dl 

ffdcOOOOOOOO 

06d2 

000000000000 

06d3 

000000000000 

06d4 

OSOOOOOOOOfO 

06d5 

000100000000 

If  [N^_„^2^1  f°'^  =max  +  10  is  not  same  (within  tolerance)  as 
+  n  =  maxrecomputeusingastill  largern.  If 

they  agree  to  go  to  0689 


Bring  up  ||  Form  abs  value 
Store  in  06D2  ||  ♦ 

Transfer  left  address  of  06C0  to  replace  right  address  of  04FF 
Go  to  069B  II  • 

Bring  up  ||  Form  abs  value 

Subtract  previous  |A^^]  1|  Form  abs  value 

Subtract  tolerance  in  06D0  II  If  A^^  is  not  precise  go  to  06CE 
If  A^^  meets  precision  go  to  064D  H  * 

Bring  up  latest  A^^  II  form  abs  value 
Store  in  06D2  ||  Go  to  069B 
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06d6 

07b3e60S8bea 

Bring  up  number  f  rom  07B3  ||  Subtract  1,000.  .  .  . 

06d7 

0bdbf406d8e2 

If  zero  go  to  06DB  11  I  f  negat  i  ve  go  to  06D8 

06dU 

0587e60000f9 

Bring  up  h  H  Absol  u  t  e  va  1  ue  o  f  h 

0bd9 

07dddc0420e2 

Is  h  larger  than  m992?  H  If  h  is  less  than  m992  go  to  0420 

Obda 

0bdee206dee2 

I  f  h  equal  s  m992  or  is  greater  thanm992  go  to  06DE 

Obdb 

0587e60000fg 

Bring  uph  11  Form  abso  1  u  t  e  va  1  ue  o  f  h 

Obdc 

07dddc0420e2 

Is  h  larger  thanm992?  H  Ifhislessthanm992  goto  0420 

Obdd 

073de2073de2 

If  h  equals  m992  or  is  greater  than  m992  go  to  073D 

Obde 

05b5eb0583ea 

Obdf 

058deb058be9 

ObeO 

07c7e807c7eb 

Obel 

07c9e807c7eb 

Obe2 

07dfe807c7eb 

0be3 

07ele807c7eb 

06e4 

07e3e807c7eb 

0be5 

07e5eb07c7eb 

Obeb 

07e7e80587fi6 

0be7 

0587eb07ece8 

Obeb 

0587eb07eee8 

Obeb 

0587eb07f0eb 

Obea 

0587eb07f2eb 

Obeb 

0383e60583eb 

Obec 

Obed 

Obee 

Obel 

ObfO 

07f4e807f4eb 

ObObebOOOOeO 

0587e607c7eb 

07f6e807c9e6 

0791e90793ec 

lines  06DE  to  0737 

compu  t  e  a  1 1 

terms  in  starting  value  of  f or  h  -•  oo 

Using  Meixner  (1)  Eq.  10  (Prolate  case) 

Obfl 

07f8e807f4e6 

0bl2 

0797eb07fae8 

06X3 

07fae707c9e9 
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06f4  0795e907c7eb 
06f5  0799ec0587ec 
0bf6  07fce807c9e6 
Otif7  079beb07ele9 
OtifS  079de90791eb 
06f9  07fee807c9e6 
Obfa  058be907f4eb 
06fb  07aleb07fae8 
Ofafc  07fae707fee9 
Obfd  U79fec07ecec 
Otife  U800e807e3e6 
Obff  07a5eb07fae8 
0700  07dfe607a7eb 
0701  07fee807c7e6 
0702  07a9eb07fee9 
0703  07fae907a3ec 
0704  07a3ec0802e8 
0705  07c7e607adeb 
0706  07fae807dfe6 
0707  07abeb07fae9 
0708  07f4eb07a3ec 
0709  0804e807c7o6 
070a  0806fib0793ec 
070b  0804ea0802e9 
070c  07eeec0808e8 
070d  07e5e607bleb 
070e  07fae807ele6 
070f  07cbeb07fee8 
0710  07c9e607cdeb 


0711 


07cfe907fee9 
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0712 

071ae907afec 

0713 

07afec0804e8 

0714 

07ele607d5eb 

0715 

07fae807c9e6 

0716 

07dleb071ee8 

0717 

07d7e907fee9 

0718 

07fRe907f4eb 

0619 

07d3ec0802e8 

071a 

07c9e6058be9 

071b 

0806eb07c5eb 

071c 

0793ec0802ea 

071d 

0804e907£0ec 

071e 

080ae807e7e6 

071f 

07d9eb07fae8 

0720 

07e3e607dbeb 

0721 

07fee807dfe6 

0722 

07b5eb0802e8 

0723 

07c7e607b7eb 

0724 

0802e9071ee9 

0725 

071ae9079fec 

0726 

079fec071ae8 

0727 

07e3e607b9eb 

0728 

071ee807dle6 

0729 

07bbeb0802e8 

072a 

07c7e607bdeb 

072b 

0802e9071ee9 

072c 

07f4et*0797ec 

072(1 

079£ec07fee8 

072e 

07b£e607d£eb 

0721 

0802e807c7e6 
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0V30 

07cleb0802e9 

U731 

0806eb07d3ec 

U732 

0802e807c7e6 

0733 

07£4eb07c3ec 

0734 

0804e80804e7 

0733 

0802c907f eea 

0736 

07£ae907f2ec 

0737 

0b0ce80865e2 

0738 

07f6e607f4e9 

0739 

07£8ea07£cea 

073a 

OfaOOeaObObea 

073b 

ObOaeaObOcea 

073c 

0obi;e80420t;2 

073d 

Oob7e605b7eb 

073e 

0b3oe80387eb 

073f 

0b3Veb05b7eb 

0740 

0b39e80o87eb 

0741 

083bc80585e6 

0742 

05b3ea03bdeb 

0743 

0583e9058be9 

0744 

083de8083deb 

0745 

0b3fe8083feb 

0746 

0b41e8083feb 

0747 

0b43e80583e6 

0748 

05b3eb0845e8 

0749 

0845Rhn847e8 

074a 

0845eb0849e8 

074b 

083de60587eb 

074c 

058dflb084be8 

074d 

083£e60845ea 
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II  Go  to  0865 


lines  0738  to  073C 

compute  starting  for  h-oo  (Prolate  case') 


II  Go  to  0420 


Lines  073D  to  078B 
com put  call 

terms  in  starting  value  of  for  h  -  oo 
Using  Meixner  (1)  Eq.  (19)  (Oblate  case) 
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07  4e 
074f 
0750 
0751 
0752 
0753 
0754 
0755 
0756 
0757 
0758 
0759 
075a 
075  b 
075c 
0754 
075e 
075  « 
0760 
0761 
0762 
0763 
0764 
0765 
0766 
0767 
0768 
0769 
076a 


058be9058dec 

084de80833ec 

0587ec083deb 

084fe80841e6 

0791eb0859e8 

0831e607ddeb 

085be8083fe6 

07c5eb058be9 

0845eb058deb 

085de8085de7 

0847e9058be9 

085be90859e9 

0799ec0835ec 

0851e80841e6 

07a5eb08Sge8 

083fe60811eb 

085be8083fe6 

080feb0815e9 

0845eb058deb 

085de80847e6 

0817eb085dea 

07abe9085be9 

0859e9083deb 

07d3ec0837ec 

0853e80843e6 

07bleb0859e8 

0841e60813eb 

085be8083fe6 

0819eb085de8 


076b 


0841e607ddeb 
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076c 

085fe8083fe6 

076d 

0801eb07c5e9 

076e 

085fe90845eb 

076f 

07ddeb085fe8 

0770 

0831e60817eb 

0771 

081de90847eb 

0772 

07c5eb0861e8 

0773 

0849e6058deb 

0774 

0863e80863e7 

0775 

0861e90851ea 

0776 

081be9085de9 

0777 

085be908S9e9 

0778 

0791ec0839ec 

0779 

0855e80843e6 

077  a 

07d9eb0859e8 

077b 

0841e6081feb 

0/7c 

085be8083fe6 

0/  /d 

0821eb083de8 

0/'7e 

0841e60825eb 

077  f 

083fe8083fe6 

0780 

0827eb0829e9 

0781 

085fe90845eb 

0782 

085fe8083Xe6 

0783 

082beb082de9 

0784 

0847eb0861e8 

0785 

0849e6082feb 

0786 

0863e80863e7 

0787 

0861e9085Xea 

0788 

0823e9085de9 

0789 

085be90859e9 
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078a  083deb0831ec 

078b  083bec086de2  H  Go  to  086D 

078c  0835e7084be9 

078d  084dea08.41ea 

lines  078C  to  0790 

078e  0851ea0853ea  compute  st  arting  value  of  for  h  -•  oo  (Oblate  case) 

0781  0855ea08S7ea 

0790  05b9e80420e2  ||  to  0420 

0791  OaOOOOOOOOOO  ) 

L - -s- 

0792  000300000000  J 

0793  080000000000 

0794  000400000000 

0795  ObOOOOOOOOOO 

0796  000400000000 

0797  080000000000 

0798  000600000000 

0799  080000000000 

079a  000700000000 

079b  OdOOOOOOOOOO 

079c  000500000000 

079d  0a8000000000 

079e  000500000000 

0791  080000000000 

07a0  OOObOOOOOOOO 

07al  OcOOOOOOOOOO 

07a2  000900000000 

07a3  080000000000 

07a4  000800000000 

07a5  084000000000 

07a6  000600000000 


07a7 


0C7400000000 
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07  aB 

OOObOOOOOOOO 

07  a9 

OafaBOOOOOOO 

07aa 

OOOdOOOOOOOO 

07  ab 

094000000000 

07  ac 

000600000000 

07  ad 

0a7000000000 

07ae 

000800000000 

07af 

080000000000 

07b0 

000900000000 

07bl 

OfcOOOOOOOOO 

07b2 

000600000000 

07b3 

000000000000 

07  b4 

000000000000 

07b5 

0£edf9000000 

07b6 

001400000000 

07  b7 

OBBdO£cOOOOO 

07bB 

001600000000 

07b9 

0b33B0000000 

07ba 

OOOdOOOOOOOO 

07bb 

0£9110000000 

07bc 

001100000000 

07bd 

091£8e000000 

07be 

001300000000 

07bX 

OblBOOOOOOOO 

07c0 

000900000000 

07cl 

0bc200000000 

07c2 

OOOf 00000000 

07c3 

080000000000 

07c4 

000500000000 

07c5 

OcOOOOOOOOOO 
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07c6  000200000000 
07c7  000000000000 
07c8  000000000000 
07cg  000000000000 
07ca  000000000000 
07cb  09a600000000 
07cc  OOOdOOOOOOOO 
07cd  0a93f0000000 
07ce  001000000000 
07cf  OafScOOOOOOO 
07d0  OOOfOOOOOOOO 
07dl  0a3c00000000 
07d2  OOObOOOOOOOO 
07d3  030000000000 
07d4  OOOaOOOOOOOO 
07d5  0e6000000000 
07d6  000700000000 
07d7  0b7c00000000 
07d8  OOOaOOOOOOOO 
07d9  083c00000000 
07da  OOOaOOOOOOOO 
07db  0X0590000000 
07dc  001000000000 
07dd  000000000000 
07de  000000000000 
07dX  000000000000 
07e0  000000000000 
07el  000000000000 
07e2  000000000000 


07e3 


000000000000 
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07  e4 

000000000000 

07e5 

000000000000 

07e6 

000000000000 

07c7 

000000000000 

07e8 

000000000000 

07e9 

00000007eal0 

07  ea 

000000000000 

07eb 

000000000000 

07  ec 

000000000000 

07ed 

000000000000 

07  ee 

000000000000 

07cf 

000000000000 

07f0 

000000000000 

07fl 

000000000000 

07f2 

000000000000 

07f3 

000000000000 

07f4 

000000000000 

07f5 

000000000000 

07f6 

000000000000 

07f7 

000000000000 

07f8 

000000000000 

07f9 

000000000000 

07fa 

000000000000 

07fb 

000000000000 

07fc 

000000000000 

07fd 

000000000000 

07fe 

000000000000 

07ff 

000000000000 

0800 

000000000000 

0801 

000000000000 
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0802 

000000000000 

0803 

000000000000 

0804 

000000000000 

0805 

000000000000 

0800 

000000000000 

0807 

000000000000 

0808 

000000000000 

0809 

000000000000 

080a 

000000000000 

080b 

000000000000 

080c 

000000000000 

080d 

000000000000 

080e 

000000080110 

0801 

0b8000000000 

0810 

000500000000 

0811 

0e4000000000 

0812 

000700000000 

0813 

OaaOOOOOOOOO 

0814 

000900000000 

0815 

OeOOOOOOOOOO 

0816 

000500000000 

0817 

OdOOOOOOOOOO 

0818 

000400000000 

0819 

OelOOOOOOOOO 

081a  000800000000 

081b  OeOOOOOOOOOO 
081c  000400000000 
081d  OeOOOOOOOOOO 
081a  000300000000 
081f  081580000000 
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0820  OOOdOOOOOOOO 

0821  0a3280000000 

0822  001100000000 

0823  0fc400000000 

0824  OOOaOOOOOOOO 

0825  OeacOOOOOOOO 

0826  OOOaOOOOOOOO 

0827  OeaGOOOOOOOO 

0828  OOOcOOOOOOOO 

0829  0c6e00000000 

082a  OOObOOOOOOOO 

082b  OeSSOOOOOOOO 

082c  000900000000 

082d  OOecOOOOOOOO 

082e  OOOaOOOOOOOO 

082f  044000000000 

0830  000600000000 

0831  080000000000 

083?.  OOOeOOOOOOOO 

0833  080000000000 

0834  000300000000 

0835  000000000000 

0836  000000000000 

0837  000000000000 

0838  000000000000 

0839  000000000000 

083a  000000000000 

083b  000000000000 

083c  000000000000 


083d 


000000000000 
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083e 

083f 

0840 

0841 

0842 

0843 

0844 

0845 

0846 

0847 

0848 

0849 

084a 

084b 

084c 

084d 

084© 

084f 

0850 

0851 

0852 

0853 

0854 

0855 

0856 

0857 

0858 

0859 

085a 


000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

oooooooooooo 

000000000000 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 


085b 


oooooooooooo 
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0S5c 

000000000000 

085d 

000000000000 

085e 

000000000000 

085f 

oooocooooooo 

0860 

000000000000 

0861 

000000000000 

0862 

000000000000 

0863 

000000000000 

0864 

000000000000 

865 

07f6c00000f9 

Bring  up  1®*  term  in  Meixner  Formula  (Prolate  case) 

II  Form  abs  value 

0866 

0835e8080ce6 

Store  in  0835  ||  Bring  up  last  term  in  Meixner  Formula 

0867 

0000f90837e8 

Form  abs  value  ||  Store  in  0837 

0868 

08t5e6086beb 

Bring  up  0835  ||  Mu  1 1 i ply  by  10”^/2 

0869 

0837ea0738f4 

Subtract  0837  ||  If  positive  go  to  0738;  If  neg.  go  to  086A 

086a 

0420e20000e0 

Go  to  0420  II  » 

086b 

0a3d70a3d70a  1 

1 

086c 

1 

tff900000000  J 

y  10-2/2  - 

086d 

0857e80835e6 

Store  in  0857  ||  Bring  up  1’^  term  in  Meixner  Formula 
(Oblate  easel 

086e 

0000f907f6e8 

Form  abs  value  ||  Store  in  07F6 

086f 

08S7e60000f9 

Bring  up  0857  (  =  last  term  in  Meixner  Formula)  ||  Form  abs  value 

0870 

07f8e807f6e6 

Store  in  07F8  ||  Bring  up  07F6 

0871 

086b<jb07f8ea 

Multiply  by  10-2/2  ||  Subtract  07F8 

0872 

078cf40420e2 

If  positive  go  to  078C;  If  negative  go  to  right  side  ||  Go  to  0420 

0873 

0025f30010,f3 

Carriage  return  ||  Space  (in  type  out) 

0874 

001013001013 

Space  J1  Space 

0875 

001013002913 

Space  II 

0876 

0610e6000bl0 

Bring  up  0610  (=tol  on  eigenvalue)  ||  Print  11 

0877 

0501e40000e0 

Enter  fixed  point  at  left  side  of  0501  |l  ♦ 

0878 

OOOOeOOOOOeO 

*  11  * 

0879 

OOOOeOOOOOeO 

•  II  • 
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087a 

0610e60791cb 

Br ing  up  Tol  of  Ei genva lue  ||  Multiply  by  5 

087b 

ObbdubOGlOcB 

Multiply  by  2  ||  Store  in  0610 

0S7c 

CGd0o80000e0 

Store  in  06D0  ||  ♦ 

087d 

OSloUOOCecOO 

Place  the  left  address  of  051E  into  the  right  address  of  069C 

087e 

0502e20000c0 

Go  to  0502  11  ♦ 

087f 

0<16(jd60CC200 

Place  the  left  address  of  046E  into  left  address  of  0692 

0880 

1002d0069200 

Decrement  left  address  of  0692  by  2 

0881 

0467c20000c0 

Go  to  0467  11  ♦ 
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